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ABSTRACT 

Context. The second parameter (the first being metallicity) defining the distribution of stars on the horizontal branch (HB) of globular clusters 
(GCs) has long been one of the major open issues in our understanding of the evolution of normal stars. Large photometric and spectroscopic 
databases are now available: they include large and homogeneous sets of colour-magnitude diagrams, cluster ages, and homogeneous data 
about chemical compositions from our FLAMES survey. 
Aims. We use these databases to re-examine this issue. 

Methods. We use the photometric data to derive median and extreme (i.e., the values including 90% of the distribution) colours and magnitudes 
of stars along the HB for about a hundred GCs. We transform these into median and extreme masses of stars on the HB, using the models 
developed by the Pisa group, and taking into account evolutionary effects. We compare these masses with those expected at the tip of the red 
giant branch (RGB) to derive the total mass lost by the stars. 

Results. We find that a simple linear dependence on metallicity of this total mass lost describes quite well the median colours of HB stars. 
Assuming this mass loss law to be universal, we find that age is the main second parameter, determining many of the most relevant features 
related to HBs. In particular, it allows us to explain the Oosterhoff dichotomy as a consequence of the peculiar age-metallicity distribution 
of GCs in our Galaxy, although both Oosterhoff groups have GCs spanning a rather large range in ages. However, at least an additional - 
third - parameter is clearly required. The most likely candidate is the He abundance, which might be different in GC stars belonging to the 
different stellar generations whose presence was previously derived from the Na-0 and Mg-Al anticorrelations. Variations in the median He 
abundance allow us to explain the extremely blue HB of GCs like NGC 6254 (=M 10) and NGC 1904 (=M 79); such variations are found to be 
(weakly) correlated with the values of the R-parameter (that is the ratio of the number of stars on the HB and on the RGB). We also show that 
suitable He abundances allow deriving ages from the HB which are consistent with those obtained from the Main Sequence. Small corrections 
to these latter ages are then proposed. We find that a very tight age-metallicity relation (with a scatter below 4%) can be obtained for GCs 
kinematically related to the disk and bulge, once these corrections are applied. Furthermore, star-to-star variations in the He content, combined 
with a small random term, explain very well the extension of the HB. There is a strong correlation between this extension and the interquartile 
of the Na-0 anticorrelation, strongly supporting the hypothesis that the third parameter for GC HBs is He. Finally, there are strong indications 
that the main driver for these variations in the He-content within GCs is the total cluster mass. There are a few GCs exhibiting exceptional 
behaviours (including NGC 104=47 Tuc and in less measure NGC 5272=M 3); however, they can be perhaps accommodated in a scenario for 
the formation of GCs that relates their origin to cooling flows generated after very large episodes of star formation, as proposed by Carretta et 
al. (2009d). 
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1. Introduction 
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mous ftp to cdsarc.u-strasbg.fr (I30.79.I28.5) or via 
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Sandage & Wallerstein (1960) noticed that the distribution with 
colour/temperature of stars on the horizontal branch (HB) of 
globular clusters (GCs) is roughly correlated with their metal 
content. A few years later, this observation was explained 
by the first successful models of HB stars describing the ef- 
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feet of metal content on the efficiency of H-shell burning in 
low mass stars where He is burning in the core (Faulkner 
1966). However, soon after this important theoretical achieve- 
ment, van den Bergh (1967) and Sandage & Wildey (1967) 
pointed out that the correlation between colour/temperature 
and metallicity had several exceptions, a difficulty that has be- 
come known as the second parameter problem. In the follow- 
ing forty years, a large number of tentative explanations for this 
discrepancy appeared in the literature, but no overall satisfac- 
tory scenario has yet been found. A proof of the large inter- 
est raised by this issue is that entering "globular cluster" and 
"second parameter" in the ADS data baseQ] resulted in 23 1 ab- 
stracts (24 since 2006) with 6031 citations on a query made 
on April 15th, 2009. Of course, this search is probably incom- 
plete, because there are many related issues, e.g., the Oosterhoff 
dichotomy in the mean periods of RR Lyrae in galactic GCs 
(Oosterhoff 1944; Sandage 1982), the UV upturn in the spectra 
of bulges and elliptical galaxies (Code 1969), the ages and the 
He abundances of GCs (Iben 1968), the mass loss law for low 
mass stars, which use different keywords. The second param- 
eter problem is certainly one of the major open issues in our 
understanding of the evolution of normal stars. For reviews of 
this topic, we refer to Moehler (2001) and Catelan (2009). 

There are various reasons why the second parameter issue 
has been insofar so difficult to solve. The most important is 
that there is most likely more than a single second parameter 
The colour of HB stars is very sensitive to several physical stel- 
lar quantities in the age and metallicity regime typical of GCs 
(see e.g.. Rood 1973; Renzini 1977; Freeman & Norris 1981; 
Fusi Pecci et al. 1993). Zinn (1980) and many authors after him 
convincingly showed that younger ages might explain the red 
colours of the HB of several of the outer halo GCs (see Dotter 
et al. 2008 for a similar line of thought). However, progress in 
the determination of the (relative) ages of GCs (e.g.. Stetson et 
al. 1996; Rosenberg et al. 1999; De AngeH et al. 2005; Marin- 
Franch et al. 2009) demonstrated that this cannot be the only 
second parameter. The same result had previously been ob- 
tained even more directly from the broad distribution in colours 
of HB stars within some individual GCs (see e.g., Ferraro et al. 
1990). 

Since a change in the mass of stars on the HB may well 
cause even large variations in their colours (Rood 1973), a 
special mass-loss law has become a popular explanation (see 
Dotter 2008 for an example of this approach). Unfortunately, 
the physics of mass loss is very poorly known at present (see 
e.g., Willson 2000; Meszeros et al. 2009; Dupree et al. 2009). 
Many different mechanisms may affect mass loss (see e.g.. Rich 
et al. 1997; Green et al. 1997; Soker & Harpaz 2000) and em- 
pirical evidence is inadequate for fully constraining them (see 
e.g., Peterson 1982; Origlia et al. 2007, 2008). Given these lim- 
itations, it is not yet possible to build an ab-initio model for 
the estimation of mass loss from GC stars. Hence we prefer to 
carefully restrict our assumptions: we looked for solutions with 
a mass loss law based on as few simple parameters as possible, 
for which we may obtain constraints from independent obser- 
vations. 



' http://adsabs.harvard.edu/ 



Additional second parameters considered included He 
abundances, the ratio of CNO to Fe abundances, stellar rota- 
tion or binarity (see e.g.. Freeman & Norris 1981), and cluster 
concentration (Fusi Pecci et al. 1993). However, all of these ex- 
planations were found to be unsatisfactory overall. In the most 
successful cases, they might explain some groups of stars with 
anomalous colours on the HB (e.g., most field O-B subdwarfs 
are binaries, see e.g., Maxted et al. 2001; Napiwotzki et al. 
2004; Han et al. 2003; however most of these in GCs seem to be 
single stars, see e.g., Moni Bidin et al. 2006). In the least suc- 
cessful cases, they are inconsistent with the data (CNO abun- 
dances: see the case of the second parameter pair NGC 362 and 
NCC 288: Shetrone & Keane 2000 and Catelan et al. 2001). 
The situations for rotation and He abundances are more com- 
plicated. For rotation, after the initial enthusiasm triggered by 
the pioneering observations of Peterson (1982), the problem 
was found to be less straightforward. More extensive data sets 
(DeHyannis et al. 1989, Behr et al. 2000a, 2000b; Behr 2003; 
Recio-Blanco et al. 2002) revealed an intricate pattern, so that 
it seems difficult to use dkect observations of rotation along 
the HB to confirm its role in the second parameter issue. While 
these observations do not rule out the possibility that rotation 
is indeed important, we cannot avoid concluding that the evi- 
dence for and against remains poor (see e.g., Sweigart 2002). 
We discuss the case of He abundances in subsequent sections. 

Over the years, an enormous wealth of observational data 
has been collected, not only in terms of the distribution of stars 
along the HB with colour, but also the chemical composition, 
the period distribution of RR Lyrae, and other properties of 
GCs. Several authors (see e.g., Fusi Pecci et al. 1993; Catelan 
et al. 2001; Recio-Blanco et al. 2006; CaiTetta et al. 2007) 
pointed out the existence of correlations between global cluster 
parameters (such as luminosity, concentration, or Galactic or- 
bit) and phenomena related to the second parameter. However, 
the mechanism connecting these global properties to the evolu- 
tion of individual stars remained elusive until a few years ago. 

A revised approach to the problem of the second param- 
eter can now be developed. It is based on what was initially 
considered to be a separate characteristic of GCs, that is the 
abundance anomalies observed for GCs in the light elements 
CNO, Na, Mg, and Al (see Kraft 1994 and Gratton et al. 2004 
for reviews of this topic). Early suggestions that there might 
be a correlation between these two sets of observations were 
made by Norris (1981), and more recently by Kraft (1994) and 
Catelan & de Freitas Pacheco (1995) and several other authors 
after them. However, the exact mechanism linking the two phe- 
nomena remained unclear. The two innovative steps taken sub- 
sequently were: 

- the recognition that typical GCs host at least two gen- 
erations of stars: this is required to explain the abun- 
dance anomalies observed for main sequence (MS) stars by 
Gratton et al. (2001) and Cohen et al. (2002). These obser- 
vations contradicted the paradigm of GCs as single stellar 
populations, and opened a new view on GC formation and 
evolution that we are now only beginning to explore (see 
Gratton et al. 2004 for early results). 
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- the understanding by D'Antona et al. (2002) that (large) 
variations in the abundance of He, which are expected to 
be correlated with the variations in Na and O and other ele- 
ments, might result in large differences in the turn-off (TO) 
masses of stars of similar age: this is because He-rich stars 
evolve faster than He-poor ones and thus, at a given age. 
He-rich stars at TO are less massive. Therefore, similar 
mass losses along the red giant branch (RGB) would lead 
to HB stars of very different masses, and hence colours. 
Not all authors estimated variations in the He abundances 
as large as those proposed by D'Antona et al. (2002; see for 
instance Marcolini et al. 2009). We note however that these 
small spreads in He abundance (which are usually justified 
on the basis of chemical evolution models) have difficulties 
in reproducing the observed spread in masses along the HB 
and the splitting of the MS of NGC 2808 (Lee et al. 2005, 
Piotto et al. 2007). 

That a combination of age and He differences may explain 
the second parameter is very attractive for several reasons: (i) 
The large variations in He and the Na-O anticorrelation are ex- 
plained by the presence of different generations of stars in GCs, 
and it is then easy to link them to general cluster properties, 
such as their mass or location in the Galaxy, which seem to 
play a role in the second parameter issue; (ii) These different 
stellar generations may well be used to explain discontinuous 
and often discrete distributions of stars along the HB, such as 
that observed e.g., in NGC 2808 (D'Antona et al. 2005); (iii) 
Very accurate photometric data detect multiple main sequences 
in some GCs (Bedin et al. 2004; Piotto et al. 2007) that can 
only be explained by assuming large variations in the He con- 
tent (Non-is et al. 2004; Piotto et al. 2005; Milone et al. 2010). 

While these observations are extremely interesting, they are 
rather limited in number: the data discussed insofar only con- 
cern a handful of massive GCs. A more comprehensive study of 
a large set of GCs, analysed in a homogeneous way, was lack- 
ing until recently. Such an analysis is now possible, thanks to 
the large databases of colour magnitude diagrams (CMD) and 
accurate ages provided by ground (Rosenberg et al. 2000a,b) 
and space observations (e.g., Piotto et al. 2002), and the ex- 
tensive data on the Na-O anticorrelation from our FLAMES 
survey (Carretta et al. 2009a,b and references therein), comple- 
mented by literature data. In this paper, we present an explo- 
ration of these unprecedented databases. In the first part of the 
paper, we consider the evidence provided by extensive photo- 
metric datasets, from both ground-based and HST observations 
for a sample of almost a hundred GCs, deriving the proper- 
ties of HBs as defined by their median values and extension, 
and examining their correlations with metallicity and age. This 
analysis produces a simple mass-loss law, that explains the me- 
dian colours of HB stars. However, as found by several au- 
thors before, an additional parameter is needed to explain the 

^ As pointed out by the referee, the knowledge that He abundance 
variations lead to variations in turn-off masses is much older (see Iben 
& Rood 1970). However, D'Antona et al. were the first to relate varia- 
tions in Na, O, and other light elements to variations in He abundances, 
and then TO masses. This created the link between abundance anoma- 
lies on the RGB and the second parameter issue. 



HB colours of GCs with an extreme blue HB (BHB) and the 
spread of colours in many other cases. In the second part of 
the paper, we consider variations in the He content as a possi- 
ble explanation of these discrepancies, we derive the He abun- 
dance variations required explaining the observed properties of 
the HBs, and discuss the implications for MS photometry. In 
the third part of the paper, we search for additional evidence 
that He is indeed the additional parameter required to explain 
the HB morphology. In particular, we explore the correlations 
with other chemical anomalies, namely the Na-O anticorrela- 
tion. For this purpose, we consider a smaller but still quite large 
sample of 24 GCs, including classical second parameter cases 
(NGC 288 and NGC 362; NGC 5272 and NGC 6205), a list 
of blue HB clusters (such as NGC 1904 or NGC 6752), and 
the very extended HB cases (such as NGC 2808, NGC 6388 
and NGC 6441). The aim of our discussion is to convince the 
reader that a combination of age and He abundance variations, 
the latter being related to multiple generations of stars within 
each GC, is a promising scenario to clarify most of the so far 
unexplained characteristics of the second parameter issue. We 
emphasise that we do consider that additional effects (e.g., bi- 
narity) may affect the colour of stars along the HB, but the roles 
played by the age and He abundance variations are probably 
dominant. 

2. Median and extreme colours of stars along the 
HB 

Several authors have suggested that age is the (main) second 
parameter determining the colour of HB stars (e.g., Zinn 1985; 
Demarque et al. 1989; Lee et al. 1994). In this section, we 
present a revised evaluation of this issue, performing an exten- 
sive comparison between observations and models. We exploit 
the large databases of globular cluster photometry available on 
the web and the latest age and metallicity estimates. 

2.1. Observational data 

The first step of our procedure was to derive the median and 
range of colours and magnitudes (hereinafter the range includ- 
ing 90% of the stars) along the HB. We derived these quantities 
from the databases of ground-based and HST (snapshot) ob- 
servations presented by Rosenberg et al. (2000a,b) and Piotto 
et al. (2002), respectively. These databases were selected be- 
cause they are publicly available and include a large number of 
clusters analysed with very homogeneous methods. The HST- 
ACS survey (Sarajedini et al. 2007) is providing new, very high 
quality data, which will largely supersede those older cata- 
logues; unfortunately, this data set is not publicly available yet. 
However, as we will see, the older photometry is of excel- 
lent quality, with only a few caveats. The ground database in- 

' An analysis of the ages of GCs and median colours of the HB 
based on these ACS data has been published after this paper was writ- 
ten (Dotter et al. 2010). A complete comparison with our results would 
require a long section. However, we note that there are very tight cor- 
relations between the V-I colours from ACS data, and the B-V and V-I 
ones of this paper. The ages by Dotter et al. also agree fairly well, at 
least in a statistical sense, with those considered in this paper. Not sur- 
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Fig. 1. Observed CMD for NGC 5904. RGB and HB stars are 
marked with red and blue points, respectively, while the green 
boxes represent the selection region of each component. 

eludes V and I CMDs and photometric tables for 52 GCs, from 
which we dropped E3, which is too scarcely populated for the 
present purposes, and u) Cen, whose spread in metal abundance 
makes the interpretation of data far more complicated than for 
the remaining clusters. The HST database contains F439W and 
F555W CMDs and photometric tables for 71 GCs. In the fol- 
lowing, we considered B and V colours obtained from the HST 
photometry through the transformations given by Piotto et al. 
(2002). There are 23 clusters in common between the two sam- 
ples, so data are available for a total of 98 GCs. 

Table[T]gives basic parameters for the clusters. Metallicities 
[Fe/H] are from our re-analysis of GC metallicity (Carretta et 
al. 2009c; note that we adopted a metallicity of [Fe/H] =-1.35 
for NGC 6256, from Stephens & Frogel 2004), while appar- 
ent distance moduli im - M)v, reddenings E{B - V), abso- 
lute magnitudes My, and HBR are from the Harris catalogue 
(1996). Relative ages are from the compilation of Carretta et 
al. (2009d), and are mainly from Marm-Franch et al. (2009) 
and De Angeli et al. (2005, this last modified to be on the same 
scale), corrected to have values consistent with the metalUcities 



prisingly, there is also agreement at least in the first major conclusion: 
age is the main second parameter. However, Dotter et al. suggest clus- 
ter concentration as third parameter, while we propose absolute mag- 
nitudes, on turn related to the He content. We note that Dotter et al. did 
not examine the extension of the HB, limiting their study to the me- 
dian colours. We argue that it is not easy to discuss this third parameter 
using this approach. Furthermore, there is some correlation between 
luminosity and concentration for GCs, which may justify their result. 
However, our analysis of the extension of the HB shows that this is 
much more strongly correlated with the cluster absolute magnitudes, 
rather than their concentration. 



Fig. 2. Comparison between the number of HB and RGB stars 
that are within ±1 mag of the mean magnitude of the HB 
V(HB). Upper and lower panels are for HST and ground-based 
data. Labels indicate the points relative to those GCs for which 
star counts are not reliable, using this diagnostics. 



by Carretta et al. (2009c). Masses at both the turn-off (Mto) 
and the tip of the RGB (Mrgb) were computed as explained in 
Sect. 12.41 IQR values are the interquartile of the Na-O, and Mg- 
Al (anti-)correlations and come from Carretta et al. (2009b), 
with a few additional data from literature (Shetrone & Keane 
2000 for NGC 362; Sneden et al. 2004 and Cohen & Melendez 
2005 for NGC 5272 and NGC 6205; Marino et al. 2008, 2009 
for NGC 6121 and NGC 6656, respectively, and Yong et al. 
2005 for NGC 6752). Finally, log T™''(HB) is the maximum 
temperature of HB stars, from Recio-Blanco et al. (2006), com- 
plemented by data for a few clusters evaluated by Carretta et 
al.(2009d). For each cluster, we first identified the region cov- 
ered by the HB, as well as that occupied by RGB stars that are 
within +1 mag of the mean magnitude of the HB V(HB), taken 
from the on line version of the Harris catalogue of GCs (Harris 
1996). Figure [U shows an example of the selection of HB and 
RGB stars. 

Completeness of the photometry can be an issue in a sta- 
tistical study such as the present one. HST-based observations 
are generally complete to magnitude V ~ 2 1 .5 (with some ex- 
ceptions on both the brighter and fainter side), which is much 
fainter than the HB level at the RR Lyrae colour in the vast ma- 
jority of GCs of this sample. We tested the completeness of the 
HB in the HST photometry by checking whether the number of 
HB stars is not unexpectedly small compared to the number of 
RGB stars that are within +1 mag of V(HB) (see Fig. |2]i. On 
average, this ratio is 0.934 ± 0.020 for GCs with HST data. The 
only cluster for which there is a clear deficiency of observed 
HB stars (a value smaller than the average one by more than 
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Fig. 4. Run of (B-V)rgb,o and (V-I)rgb,o with metallicity. 
Open and filled symbols represent clusters with E(B-V) larger 
and smaller than 0.25, respectively, while the solid curves show 
the relations given by Equations ([T) and (|2]i. A few points rela- 
tive to interesting clusters are marked. 

twice the standard deviation, as estimated by Poisson statis- 
tics) is the farthest one (NGC 2419: Harris et al. 1996), for 
which this ratio is 0.56 + 0.09. The next smallest value is for 
NGC 6624, for which the value is 0.70 + 0.14, which is only 
1.6 standard deviation below the average value. However, the 
ground-based observations are shallower, with typical limiting 
magnitudes in the range 19-20. For this reason, extreme HB 
stars are sometimes missing from these data; we then generally 
give preference to the HST snapshot data in our discussion. 
However, the ground-based data provide very useful informa- 
tion about a number of nearby GCs missing in the Piotto et al. 
(2002) databasfl 

Particular care was devoted to the separation between the 
HB and the RGB, which is difficult for very metal-rich GCs 
when differential reddening is large. We also tried to mini- 
mize the impact of contamination by field stars, which may 
significantly affect the determination of the colour extremes, 
but only marginally affect the determination of median colours 
and magnitudes, which are robust estimators. Unfortunately, 
while very extensive and homogeneous, the databases that we 
used include neither comparison fields allowing a (statistical) 
subtraction of field contaminants, nor a full membership study. 
Our procedure was to first estimate those field stars expected 

The databases we considered actually also include entries for vari- 
able stars (on the HB mainly RR Lyrae). These entries are very inaccu- 
rate, since usually data at very few epochs are available, and photom- 
etry in different bands is not based on data acquired simultaneously. 
However, while inaccurate, these data still have some statistical mean- 
ing; for this reason, they are considered throughout our analysis. 



to fall within the region of the CMD that we identified with 
the HB (Col. 3 of Table |2l) and the RGB, based on the galac- 
tic model TRILEGAL (Girardi et al. 2005; Vanhollebeke et al. 
2009). We then subtracted from the CMD the star closest to 
each field star (weighting differences in colours five times more 
than differences in magnitudes); the numbers of HB and RGB 
stars given in this paper are those after this subtraction. This 
procedure introduces some uncertainty in our results, in partic- 
ular in cases of strongly contaminated fields. However, eye in- 
spection of the colour-magnitude diagrams obtained using this 
procedure shows that it worked satisfactorily. Figure[3]shows a 
couple of examples of this decontamination. 

We typically identified a few hundred HB stars, and a simi- 
lar number of RGB stars. However, in a few cases we identified 
only a few tens of stars belonging to the two sequences. In these 
cases, larger uncertainties are associated with the extreme val- 
ues, while the median is still quite robust (see however Sect.|6]l. 

We fitted each observed HB with a polynomial (using the 
colour as an independent variable), whose degree varied (be- 
tween 1 and 4) from cluster-to-cluster, depending on the exten- 
sion of the HB. For each star, we then replaced the observed 
colour and magnitude with that corresponding to the closest 
position along the polynomial. The metrics that we consid- 
ered weight differences in colours five times more than those 
in magnitudes. This procedure reduces the impact of outliers. 
We then ordered the stars in terms of either increasing colour or 
magnitude, and determined the median as well as the colours 
and magnitudes that include the central 90% of the distribution. 

In the case of the RGB, we determined the average colour 
of the sequence at the magnitude of the HB by fitting a straight 
line through the points corresponding to the RGB-selected 
stars; in this case, we used the magnitude as an independent 
variable. 

Values of observed colours and magnitudes are given in 
Tables |2] and [3] for HST and ground-based data, respectively. 
Error bars are only those from statistics. They do not include 
systematic errors, which are due to incompleteness at faint 
magnitudes (most relevant for ground-based observations), in- 
correct separation of the different sequences (sometimes possi- 
ble for most metal-rich clusters with differential reddening), or 
uncertainties in decontamination by field stars. 

2.2. From observed to intrinsic colours 

Observed colours and magnitudes differ from intrinsic ones for 
mainly three reasons: (i) errors in photometric calibration and 
transformation to the standard system (colour and magnitude); 
(ii) distance to the stars; and (iii) interstellar absorption and red- 
dening. To reduce the impact of these uncertainties, we used the 
difference between observed and predicted colours of the RGB 
at the HB level ((B - y)RGB,o and (V - /)rgb,o for HST and 
ground-based observations, respectively) to estimate the offset 
in colours appropriate for each cluster, and then corrected the 
median and the extreme colours of the HB for these offsets. To 
a first approximation, errors in reddening and photometric cal- 
ibration should cancel out when determining these quantities. 
However, (B - y)RGB,o and (V - /)rgb,o both depend on metal- 
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Fig. 3. Two examples of the application of the statistical field subtraction used in this paper (top: NGC 3201 ; bottom NGC 6779); 
both cases are based on ground-based photometry. The left panels give the original photometry; the central panels give the CMD 
of field stars appropriate for the direction and reddening of the GCs, from the TRILEGAL model of the Milky Way (Girardi et 
al. 2005; Vanhollebeke et al. 2009); the right panels are the field-subtracted HB (red points) and RGB (blue points). The RGB 
stars are only those within +1 mag of the HB level, according to Harris (1996). 



licity. To estimate this dependence, we fitted with polynomials 
the relations obtained for clusters with only moderate redden- 
ing (E{B - V) < 0.25), using reddening from Harris (1996) 
and metallicities from Carretta et al. (2009c). We obtained the 
relations 

(B - y)RGB,o = 1 -235 + 0.4397 [Fe/H] + 0.0927 [Fe/H]^ (1) 
and 

(V - /)rgb,() = 1-154 + 0.2066 [Fe/H] + 0.0400 [Fe/H]^ (2) 

for HST (33 clusters, r.m.s. =0.026 mag) and ground-based data 
(29 clusters, r.m.s. =0.039 mag), respectively (see Fig.|4]i. 



Inspection of the upper panel of Fig. |4] reveals that there 
is a systematic trend in which (B - y)RGB.{) values for highly 
reddened metal-rich GCs are systematically below the calibra- 
tion curve, when using the HST photometry. This discrepancy 
might have either of three causes (or a combination of them): 
(i) metallicity or (ii) reddening for these clusters are overes- 
timated; or (iii) the colour transformation at very red colours 
((B - V) > 1.5) is incorrect, producing too red colours for 
these stars. Errors in reddening have no impact in our analy- 
sis, since the procedure we adopted should cancel their effect. 
The impacts of errors in the metallicity scale and colour trans- 
formations are more difficult to trace throughout the analysis. 
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Table 1. Basic parameters for globular clusters considered in this paper (Mto and Mrgb are in solar masses). The complete table 
is available only in electronic form. 



Cluster 


Other 


[Fe/H] 


(m - M)v 


E{B - V) 


Mv 


HBR 


Age 


Mto 


Mrgb 


IQR[Na/0] 


IQR[A1/Mgj 


logTX(HB) 






dex 


mag 


mag 


mag 






Mo 


Mo 


dex 


dex 


K 


NGC 104 


47Tuc 


-0.76 


13.37 


0.04 


-9.42 


-0.99 


0.95 


0.862 


0.909 


0.472 


0.091 


3.756 


NGC 288 




-1.32 


14.83 


0.03 


-6.74 


0.98 


0.90 


0.827 


0.869 


0.776 


0.059 


4.221 


NGC 362 




-1.30 


14.81 


0.03 


-8.41 


-0.87 


0.80 


0.858 


0.900 


0.670 




4.079 


IC 1257 




-1.73 


19.25 


0.73 


-6.15 


-0.71 














NGC 1261 




-1.27 


16.10 


0.01 


-7.81 


1.00 


0.79 


0.847 


0.887 






4.079 



Table 2. Photometric data from the HST-snapshot (Piotto et al. 2002). The complete table is available only in electronic form. 



Cluster 


n(HB) 


Npicld 




(F439 - F555)„ed 


F555 


(F439 - F555),„„ 


(F439 - F555)„,„ 


n(RGB) 


(F439 - F555)rgb 








mag 


mag 


mag 


mag 


mag 




mag 


NGC 104 


363 





14.062+0.002 


0.815+0.005 


14.047+0.003 


0.719+0.005 


0.896+0.081 


432 


1.067+0.004 


NGC 362 


274 





15.517+0.003 


0.563+0.005 


15.900+0.168 


-0.003+0.036 


0.691+0.051 


282 


0.865+0.004 


IC 1257 


34 


3 


20.440+0.063 


0.722+0.032 


20.645+0.031 


0.599+0.024 


0.889+0.028 


38 


1.522+0.013 


NGC 1261 


135 





16.821+0.009 


0.584+0.011 


17.324+0.129 


-0.047+0.038 


0.697+0.022 


124 


0.849+0.006 


NGC 1851 


307 





16.205+0.017 


0.572+0.026 


16.908+0.094 


-0.057+0.010 


0.705+0.088 


310 


0.919+0.004 



Table 3. Photometric data from ground-based data (Rosenberg et al. 2000a,b). The complete table is available only in electronic 
form. 



Cluster 


n(HB) 


Npicw 


Vmed 


{V - /)n,ed 


V ■ 

' min 


(V - Dm,. 


(V - /)n,ax 


n(RGB) 


{V - /)rgb 








mag 


mag 


mag 


mag 


mag 




mag 


NGC 104 


192 





13.970+0.002 


0.892+0.006 


13.996+0.003 


0.817+0.009 


0.955+0.004 


225 


1.054+0.005 


NGC 288 


83 





16.083+0.114 


0.022+0.022 


17.110+0.129 


-0.136+0.016 


0.228+0.531 


101 


0.987+0.006 


NGC 362 


56 


1 


15.371+0.010 


0.779+0.015 


15.523+0.003 


0.135+0.020 


0.838+0.034 


69 


0.973+0.007 


NGC 1261 


221 


1 


16.653+0.003 


0.736+0.005 


16.942+0.272 


0.026+0.095 


0.826+0.021 


239 


0.926+0.004 


NGC 1851 


33 


3 


16.177+0.020 


0.800+0.037 


17.086+0.143 


-0.028+0.014 


0.855+0.009 


35 


1.004+0.009 



However, in general they would lead to a small underestimate 
of the masses of HB stars in these clusters; this would reduce 
the sensitivity to metallicity of the mass-loss law we derive 
later, but should have no other major impact on our analysis. 

We estimated intrinsic colours by adding the offsets from 
the mean relations (for (B - y)RGB,o and {V - /)rgb,o) to the 
observed colours. Median and extreme of the colours are plot- 
ted against metallicity in Figs. |5] and |6] for HST and ground- 
based data, respectively. Different symbols are used for clus- 
ters of different relative ages (see CaiTetta et al. 2009d). In 
particular, we note that old clusters define a tight, unique (al- 
though quite complex) relation between metallicity and median 
colour of the HB. This indicates that the median location of the 
HB is uniquely determined by age and metallicity. The rela- 
tion between colours and masses of HB stars is however not 
linear, with the well-known, very strong sensitivity of colours 
on metallicities for the metallicity range -1.6 <[Fe/H]< -1.1. 
This is more clearly explained in the next subsection. Younger 
clusters typically have redder median colours than older ones 
at a given metallicity, as expected. Finally, clusters with no age 
determination roughly occupy the same region identified by the 
other clusters, but they seem to have a larger scatter. This is not 
unexpected, since data for these GCs are of poorer quality, and 
this is the reason why they lack an age estimate. We note also 
that much larger scatters are obtained for both minimum and 
maximum colours. We return to this point in the following. 



In the case of magnitudes, we simply corrected the ob- 
served values for the apparent distance moduli given in Harris 
(1996). 

We note that the plot of the median colour with metal- 
licity is conceptually similar to a Lee-Zinn diagram (Lee & 
Zinn 1990), since the median colours are very well correlated 
with the HBR parameter used in this diagram (see Fig. |7J for 
a definition of HBR, see Zinn 1985). The only exception to 
this close coiTelation is NGC 2808, when using HST (but not 
ground-based) data. This is due to the multi-modal distribu- 
tion of stars along the HB of this cluster, and to the incom- 
pleteness of the ground-based photometry, which was used to 
derive the value of the HBR parameter. On the other hand, min- 
imum colours and maximum absolute magnitudes along the 
HB are obviously well correlated with the maximum temper- 
atures logT™''(HB) determined by Recio-Blanco et al. (2006, 
see Fig. [S), which were derived from the same HST photome- 
try. The few discrepant cases are either GCs with poor statistics 
or CMDs that are heavily contaminated with field stars, and for 
which different assumptions about the impact of contaminants 
were adopted here with respect to Recio-Blanco et al. (2006). 

Values of intrinsic colours and magnitudes are given in 
Tables |4]and|5]for HST and ground-based data, respectively. 
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Table 4. Intrinsic HB photometric parameters for GCs with HST data. The complete table is available only in electronic form. 



Cluster 


(B - V)o.RGB 


(B - 


(5- V)„,ed 


(S - vu. 


Mvmax 


Mymed 




mag 


mag 


mag 


mag 


mag 


mag 


NGC 104 


0.968 


0.652 


0.736 


0.807 


0.651 


0.665 


NGC 362 


0.803 


-0.017 


0.547 


0.665 


1.090 


0.684 


IC 1257 


0.649 


-0.029 


0.083 


0.231 


1.371 


1.164 


NGC 1261 


0.809 


-0.044 


0.587 


0.690 


1.227 


0.697 


NGC 1851 


0.860 


-0.096 


0.534 


0.656 


1.442 


0.712 



Table 5. Intrinsic HB photometric parameters for GCs with ground-based data. The complete table is available only in electronic 
form. 



Cluster 


(V - /)o.RGB 


iV - /)„„„ 


{V - /)„,ed 




Mymax 


M^med 




mag 


mag 


mag 


mag 


mag 


mag 


NGC 104 


1.003 


0.783 


0.858 


0.921 


0.626 


0.600 


NGC 288 


0.949 


-0.172 


-0.014 


0.192 


2.280 


1.253 


NGC 362 


0.935 


0.115 


0.759 


0.818 


0.713 


0.561 


NGC 1261 


0.913 


0.056 


0.766 


0.856 


0.842 


0.553 


NGC 1851 


0.978 


-0.066 


0.762 


0.817 


1.616 


0.707 
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Fig. 5. Median and extreme colours as a function of metallic- 
ity for clusters observed with HST. Filled (red) triangles and 
(blue) dots indicate GCs older and younger than a relative age 
parameter of 0.92, respectively; clusters for which an age es- 
timate is not available are marked with (black) empty squares. 
Points representing a few interesting clusters are marked. 



2.3. From intrinsic colours to masses 

Our next step was to transform observed colours into masses 
of stars along the HB. The masses were obtained by comparing 
the observed colours and magnitudes with predictions from HB 
evolutionary models. For this purpose, we used the database of 



-2.5 -2 -1.5 -1 -0.5 

[Fe/H] 

Fig. 6. Same as Fig.|5]but for ground-based sample. 



models computed by the Pisa evolutionary group (Castellani 
et al. 2003, Cariulo et al. 200411 which are particularly useful 
here since they provide a grid of HB evolutionary sequences for 
different masses and metallicities. We first derived the masses 
appropriate for ZAHB stars of the same colours/magnitudes of 
the observed loci on the HB, and then applied corrections ap- 
propriate to taking into account the evolution of the stars away 
from the ZAHB. 



5 http://astro.df.unipi.it/SAA/PEL/ZO.html 
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Fig. 7. Median colours versus HBR ratio for HST (upper panel) 
and ground-based (lower panel) clusters. Different symbols de- 
note different ages (see Fig.|5]l. 



-1-0.1425 (B-V)- 0.6560 (B - vf + 

0.6277 (B - V)^ M/Mq, (3) 

and 

M(ZAHB) = 0.5271 - 0.0629 [Fe/H] - 0.0851 [Fe/H](y - /) 
-hO.1399 (V - /) - 0.6490 (V - if 
-h0.6134(y-/)3 M/Mq, (4) 

which are valid in the range 0.54 < Mzahb < 0.72 M/Mg and 
-2.5 <[Fe/H]< -0.6 (some extrapolation is required for the 
most metal-rich GCs). The similar equation for the V magni- 
tude is 

M(ZAHB) = 0.5544 - 0.07286 [Fe/H] + 0.005228 [Fe/H] 

-0.04555 My H- 0.010522 M/M^, (5) 

which is vaHd for My > 0.8 and -2.5 <[Fe/H]< -0.6. 

The evolutionary corrections were obtained by deriving 
masses from median and extreme colours of synthetic HB dia- 
grams, obtained by taking into account evolution, with respect 
to the values appropriate for the same set of stars when on the 
ZAHB. We obtained the correcting formulae 

M™i„(ev) = M„in(ZAHB) - (0.0331 - 0.0140 [Fe/H] 

-0.099 M^n(ZAHB)) M/M^, (6) 



> 



• • _ ,NGC638a 



H — H 



• ■ • 



NGC6235, 



Mined(ev) 



and 



log (HB) (K) 

Fig. 8. Run of minimum colours and maximum magnitudes 
with the maximum temperatures logT™''(HB) retrieved by 
Recio-Blanco et al. (2006). 



For the first step, we derived transformations from the evo- 
lutionary tracks given by 



M(ZAHB) = 0.5254 - 0.0650 [Fe/H] - 0.1181 [Fe/H](B - V) 



M„ed(ZAHB) - (0.0246 - 0.0136 [Fe/H] 
-0.0733 M„ed(ZAHB)) M/Mg, 



(7) 



Mn,ax(ev) = M„,ax(ZAHB)- [0.0356 -H 0.1972 [Fe/H] 

H-0.1020 [Fe/nf - M„ax(ZAHB)(0. 18936 [Fe/H]) 



H-0.10727[Fe/H]2] M/Mg 



(8) 



for minimum, median, and maximum mass, respectively. 

In practice, we estimated masses from: (i) colours when- 
ever these were redder than (B - V)o > 0.1 (or (V - /)o > 0.2), 
(ii) magnitudes whenever colours were bluer than (B - V)o = 
(V - /)o < -0.1, and (iii) by a weighted average of the values 
obtained from magnitudes and colours when they were in an 
intermediate range. 

Not unexpectedly, there is a good correlation between the 
differences M^ax - M^i„ and M,ned - Mmin (see Fig.|9]l, the re- 
gression line (for HST data) being M^ax - M„iin = (1.184 + 
0.097) (Mn,ed - M„i„) + (0.033 + 0.026 Mq with a linear re- 
gression coefficient of r=0.82 (74 clusters). We note, how- 
ever, that the estimate of the maximum mass is of much 
larger uncertainty than that of the median mass, because of the 
large evolutionary effects and the impact of smaller statistics. 
Furthermore, the median mass is on average much closer to the 
maximum than to the minimum mass, because of the difference 
between these two values is mainly due to a nearly constant off- 
set of about 0.03-0.04 Mq. In only a few extreme cases (e.g., 
NGC 2808) is the difference between M^ed and M^ax suffi- 
ciently large to really affect the arguments made in the remain- 
der of this paper Hereinafter we thus do not refer to M^ax in 
our discussion. 
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Table 7. Masses on HB for GCs with ground-based data. The 
complete table is available only in electronic form 



Cluster 




M„ed 




AM„ax 


AM,^ed 




Mo 


Mo 


Mo 


Mo 


Mo 


NGC 104 


0.650±0.003 


0.674±0.003 


0.691 


0.265 


0.242 


NGC 288 


0.570±0.001 


0.603±0.004 


0.645 


0.302 


0.269 


NGC 362 


0.638±0.004 


0.702±0.006 


0.717 


0.262 


0.197 


NGC 1261 


0.625±0.004 


0.701±0.005 


0.730 


0.279 


0.204 


NGC 1851 


0.586±0.003 


0.688±0.009 


0.702 


0.320 


0.219 



Fig. 9. Differences between maximum and minimum masses 
along the HB, compared to differences between median and 
minimum masses, for clusters observed with HST. Filled (red) 
triangles and (blue) dots indicate GCs older and younger than a 
relative age parameter of 0.92, respectively; clusters for which 
an age estimate is not available are marked with (black) empty 
squares. 

Table 6. Masses on HB for GCs with HST data. The complete 
table is available only in electronic form 



Cluster 


M ■ 


Mn,ed 


M 


AM„,,„ 


AM,„ed 




Mo 


Mo 


Mo 


Mo 


Mo 


NGC 104 


0.629±0.002 


0.648±0.001 


0.666 


0.286 


0.268 


NGC 362 


0.602±0.002 


0.680±0.011 


0.704 


0.297 


0.219 


IC 1257 


0.615±0.011 


0.655+0.001 


0.685 






NGC 1261 


0.593+0.001 


0.683+0.007 


0.707 


0.311 


0.222 


NGC 1851 


0.579+0.003 


0.664+0.005 


0.687 


0.328 


0.243 



We note that the lowest mass HB models considered here 
have a mass of 0.52 Mq, the corresponding absolute magnitude 
predicted for these stars being My < 4. However, there are a 
few GCs for which the maximum absolute magnitude of HB 
stars is My > 4. In these cases, the masses given by our equa- 
tions may be lower than the core masses at the He-flash for 
these models (Cassisi et al. 1998). While this might be indica- 
tive of either a different origin or a late He-flash for these stars 
(see e.g., D'Cruz et al. 1996), we note that in these cases our 
mass estimates are inferred from an extrapolation outside the 
range covered by the stellar models, and may therefore have a 
larger uncertainty. The value of the minimum mass along the 
HB for these GCs is certainly small, although probably not as 
small as indicated by our formulae. 

Statistical errors in the mass estimates can be obtained by 
combining statistical errors in the colours with the errors in the 



metallicities. Errors in colours cause an uncertainty of about 
0.007 Mq in the masses, while masses change by -0.005 Mq if 
metallicity is increased by -i-O.l dex in [Fe/H]. 

The comparison between HST and ground-based results is 
as follows: 

- Minimum mass: ground-HST=0.025 + 0.005Mq (rm.s 
0.024 Mq, 27 clusters), where the scatter is mainly caused 
by GCs (NGC 2808, NGC 5986, and NGC 6266=M 62) 
with (m-M)v > 15.5 and with BHB tails that are fainter 
than the limiting magnitude of the ground-based photome- 
try. When these three clusters are eliminated, the offset is 
0.019 + 0.003 Mq (r.m.s 0.016 Mq, 24 clusters), 

- Median mass: ground-HST=0.019 + 0.003 Mq 
(rm.s=0.019 Mq, 28 clusters). The only clear out- 
Her is NGC 2808. The off^set is 0.016 + 0.003 Mq 
(rm.s=0.015 Mq, 27 clusters) if this cluster is eliminated. 
This confirms that the median mass is far less sensitive 
to outliers and photometric errors than the minimum 
mass. The case of NGC 2808 is very peculiar, because 
the distribution along the HB is very discontinuous and 
clumpy. If the extreme blue tail were neglected (as in the 
case of the ground based photometry), the median would 
jump from the blue HB to the red HB, with a change in 
mass of about 0.04 Mq. 

Figure [TO] illustrates the run of median and minimum 
masses with metallicity [Fe/H], from HST and ground-based 
data, respectively. As for colours, we indicate with different 
symbols clusters of different relative ages. If we consider only 
the old clusters, the plots with median masses show a quite 
simple (although not linear) relation with metallicity, with a 
very small scatter about the mean line. Younger clusters scatter 
above the relation defined by the old ones, as expected. Again, 
clusters with no age estimates occupy a similar locus, but with 
a considerable scatter 

As noted in the case of colours, when we also consider min- 
imum masses the scatter is much larger than obtained for the 
median. We underline that in spite of some uncertainties in our 
mass derivations this scatter is real, corresponding to a large 
variation in the appearance of the HBs: we note for instance 
the low minimum masses that are obtained for the GCs with 
very extended blue tails, such as NGC 2808 or NGC 7078. We 
return to this point in Sect. |6] 

Values for the minimum, median, and maximum masses 
obtained following the previous procedure are given in Tables|6] 
and|7] for HST and ground-based data, respectively. 
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Fig. 10. Median and minimum masses vs. [Fe/H] for HST and ground based samples. Symbol legend is given in Fig. |5] A few 
points relative to interesting clusters are marked. 



2.4. Mass-loss law 

The tight relation between median HB mass and metallicity for 
old GCs obtained in the previous subsection suggests that the 
median location of stars along the HB of a GC can be deter- 
mined using a combination of metallicity, age, and a relatively 
simple and uniform mass-loss law. To show that this is indeed 
the case, we first estimated the original masses of stars cur- 
rently on the HB. To determine this value, we first consider 
masses at TO (Mjo), which can be obtained from their age, 
chemical composition, and He content, by using appropriate 
models. In practice, we used the relation 

Mto = 10""-^^'"°*5^°'=(0.9636-H 0.171 [Fe/H] 

H-0.04073 [Fe/H]2) Mq (9) 

that closely follows the isochrones of the Pisa group when 
Y=0.25, over the range of parameters (age, metallicity) appro- 
priate for GC stars. We note that Age in this formula is the rel- 
ative age parameter of Carretta et al. (2009d), and we assumed 
that Age=l corresponds to an isochrone age of 12.5 Gyr We 
then corrected this value to those appropriate for stars at the 
tip of the RGB (neglecting mass loss, hence assuming the mass 
of the stars to be identical to those at the beginning of their 
evolution), using the formula 

Mrgb = Mto + 0.056 + 0.01 17 [Fe/H] Mq, (10) 



which again fits the Pisa group isochrones. 

We then assume that the original mass of stars at the tip 
of the RGB is equal to that of stars currently on the HB (on 
average, it should be slightly lower, but this difference can 
be neglected in our discussion). The median mass lost by the 
stars before reaching the HB can then be obtained by com- 
paring their current mass Mmed (as determined in the previous 
paragraph) with Mrgb- The precise evolutionary phase (before 
ZAHB) at which mass loss occurs is not important in our ap- 
proach. It would of course play an important role in the com- 
parison with mass loss laws. The values of the median mass 
loss for stars in each cluster are plotted against metallicity in 
Fig.[TT]for clusters with HST (upper panel) and ground-based 
data (bottom panel). As usual, we used different symbols for 
old and young clusters (age is required when deriving this es- 
timate of the mass loss, so that clusters with no age estimates 
cannot be plotted here). Mass losses obtained from the two sets 
of data are quite similar, with a small offset of about 0.020 Mq 
between them. The best-fit straight line through the HST data 
is 

M,os, = (0.359 + 0.024) + (0.0942 + 0.0066) [Fe/H] Mg, (1 1) 

which has a highly significant correlation coefficient of r=0.91 
over 45 clusters. The scatter of points relative to each cluster 
around this mean line is small (rm.s. of 0.021 and 0.027 Mq 
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Fig. 11. Median mass loss as a function of metallicity for GCs older (triangles) and younger (dots) than 0.92. The solid line is 
the best-fit relation to the HST data, while the dashed one is just shifted by 0.02Mo, reproducing rather well the ground-based 
data-set. A few points relative to interesting clusters are marked. 



for old/young clusters, respectively) from HST data. The scat- 
ter is very similar (0.022 and 0.025 M© for old/young clusters, 
respectively) for ground-based data. 

At first look, this small scatter appears to be compati- 
ble with errors in metallicity and ages that are used in this 
derivation. For instance, in the case of the metal-rich cluster 
NGC 6441, changing the assumed values of the metallicity 
[Fe/H] and age parameter by 0.10 would cause variations in 
the values estimated for the mass loss by 0.014 and 0.031 M© 
(larger mass losses are obtained if metallicity and ages are in- 
creased). Sensitivities to these parameters are smaller for more 
metal-poor clusters; e.g., in the case of NGC 7078, similar 
changes in the assumed values of metallicity and age would 
cause variations in the mass loss of 0.007 and 0.022 M©. In con- 
clusion, a simple unique mass-loss law, combined with metal- 
licity and age, seems able to define quite well the median mass 
and colours of HB stars for the whole sample of clusters con- 
sidered in this analysis. However, as we will show in Sect. |4] 
the deviations of individual points from the mean relation are 
caused mostly by true physical effects. 

We may compare this estimate for the total mass lost along 
the RGB with predictions by mass-loss laws in the literature. A 
compilation of these predictions is given in Table 1 and Fig. 4 
of Catelan (2009). The present estimate of the mass lost is gen- 



erally much higher than the values given by these predictions; 
only a small fraction (0.002 M©) of this difference can be at- 
tributed to the identification of the mass of HB stars with that 
of stars at the tip of the RGB. The observed run compares quite 
well with a simple or modified Reimers (1975a,b) laws, or the 
law by VandenBerg et al. (2000), provided that the efficiency 
parameter is roughly doubled. On the other hand, other mass- 
loss laws (Mullan 1978; Goldberg 1979; Judge & Stencel 1991) 
provide a total mass lost along the RGB that changes too fast 
with metallicity to reproduce current results. 

Origlia et al. (2007) proposed a mass-loss law based on 
mid-lR Spitzer observations of red giants in globular cluster. 
According to this mass-loss law, GCs star should lose about 
0.23 Mq while ascending the RGB, nearly independent of their 
metallicity. While this average value roughly corresponds to 
what is needed to explain the HB of GCs of intermediate metal- 
licity, our analysis indicates that the total mass lost along the 
RGB should depend on metallicity. 

We note that the assumption of a particular age-metallicity 
relation (the one by Marin-Franch et al.) is implicit in our anal- 
ysis. Had we adopted a different age-metallicity relation, the 
metallicity-mass-lost relation would be different. For instance, 
if the (relative) ages of disk-inner halo clusters were inde- 
pendent of metallicity (rather than with a slope of ~ -0. 1 as 
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Fig. 12. RR Lyrae richness parameter, S(RR), as a function of 
the RR Lyrae fraction for both HST (upper panel) and ground- 
based sample (lower panel). 



adopted throughout this paper), the slope of the mass-loss law 
with metallicity would be reduced by ~ 0.02 M^jdex, from 
~ 0.09 M^/dex to ~ 0.07 M^ldex. This would not modify sig- 
nificantly any of the main results of the paper. 

3. RR Lyrae 

RR Lyrae variables are HB stars within the classical instability 
strip. Properties of the pulsating stars may be used to constrain 
several basic properties, including mass, radius, and chemical 
composition. An in-depth star-by-star comparison, while very 
illuminating, requires a considerable effort, because individual 
HB stars may be within the instability strip at various phases 
of their HB evolution. This comparison then requires a synthe- 
sis of populations, such as those performed by Marconi et al. 
(2003), Di Criscienzo et al. (2004), Caloi & D' Antona (2007, 
2008), and D' Antona & Caloi (2008). We defer this analysis 
to forthcoming papers; here we briefly examine only a few fea- 
tures. 

3.1. Data 

Because of the uncertainty related to evolution, we focus only 
on those clusters where the RR Lyrae population is dominated 
by stars still close to the ZAHB in the CMD. Since most of the 
lifetime of HB stars is spent close to their ZAHB location, we 
may be reasonably sure that this occurs for GCs that contain 
a large population of RR Lyrae. To identify these clusters, we 
consider the RR Lyrae richness parameter S(RR) as tabulated in 
the Harris (1996) catalogue, which is the number of RR Lyrae 



per unit cluster luminosity. llnFig.m we compared this pa- 
rameter with the number of stars in the instability strip that 
we may obtain from the photometric catalogues considered in 
Sect. |2] We derived the number of RR Lyrae in each cluster by 
assuming that all HB stars whose dereddened (B - V\) colour 
is 0.15 <{B- V\) < 0.45 + 0.024 ([Fe/H] + 1.5) are RR Lyrae 
(see Marconi et al. 2003 and Di Criscienzo et al. 2004). 

In Table[8] we report for each cluster the S(RR) values (Col. 
2) as well as the number and the fraction of RR Lyrae stars 
from both HST and ground-based samples (Cols. 3, 4, 5, 6). In 
Col. 7, we list the f(RR) estimates derived whenever possible 
from S(RR) values, or otherwise from f(RR), while the mean 
period (Pab) and the corresponding number of stars (n(RRab)) 
are reported in Col. 8 and 9, respectively. These last two values 
are from Clement et al. (2001). 

We note that the photometry considered throughout this 
paper is based on only a few, not even simultaneously ac- 
quired images in each filter. Hence, individual RR Lyrae will 
have nominal colours and magnitudes that can differ substan- 
tially from their mean values, and may even locate them out 
of our definitions of both the instability strip and the HB lo- 
cus. However, we might expect that, provided that photomet- 
ric errors are not too large, there should be a correlation be- 
tween the number of stars that according to the photometric 
catalogues are within the instability strip, and the true number 
of RR Lyrae. This appears to be the case: the correlation be- 
tween the fraction of HB stars that are nominally within the 
instability strip, and the RR Lyrae richness parameter listed by 
Harris, is quite good for HST data which have smaller photo- 
metric errors. It is somewhat poorer, but still fairly good for 
ground-based data. This correlation allows us to (i) calibrate 
the richness parameter in terms of the fraction of HB stars that 
are actually RR Lyrae f(RR) (we found that this is represented 
by the mean relation f{RR) = 0.01 -t- 0.00458 S{RR)); and (ii) 
complement the RR Lyrae richness parameter with data ob- 
tained from the photometry of those few clusters where this 
datum was missing. 

To estimate the separation between clusters whose RR 
Lyrae population is dominated by ZAHB objects rather than 
evolved ones, we considered in more detail a few interesting 
cases. We found that using parameters that provide an accu- 
rate distribution of stars along the HB, the fraction of RR Lyrae 
stars that were within the instability strip also when they were 
on the ZAHB is 0.67 for NGC 4590 (f(RR)=0.23), 0.45 for 
NGC 7078 (f(RR)=0.10), and close to zero for NGC 4833 
(f(RR)=0.04). We note however that the values of f(RR) are 
quite uncertain, and that synthetic HBs need several param- 
eters to reproduce various aspects of the colour distribution. 
The fractions cited above are then uncertain and should only 
be taken as indicative of a decreasing ratio of stars close to 
the ZAHB to the total of RR Lyrae variables with decreas- 
ing specific frequency; and that a value in the range of f(RR)~ 
0.04 -0.1 roughly corresponds to the transition between a pop- 



^ As pointed out by the referee, the value of S(RR) listed in the 
Harris catalogue are underestimated for some GCs. A couple of ex- 
amples ai-e NGC 2808 (see Corwin et al. 2004) and M 62 (Contreras 
et al. 2005). 
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Table 8. Properties of the RR Lyrae variables for each cluster The complete table is available only in electronic form. 



Cluster 


S(RR) 


n(RR)HST 


f(RR)HST 


n(RR)GB 


f(RR)GB 


f(RR) 


<Pab > 

days 


n(P„b) 


NGC 104 


0.2 


2 


0.01 


1 


0.01 


0.01 


0.737 


1 


NGC 288 


4.0 






3 


0.04 


0.03 


0.678 


1 


NGC 362 


3.0 


35 


0.13 


4 


0.07 


0.02 


0.542 


7 


IC 1257 







0.00 






0.00 






NGC 1261 


13.5 


8 


0.06 


21 


0.10 


0.07 


0.555 


13 



Table 9. Comparison of median and maximum masses with the 
masses for RR Lyrae variables derived by Di Criscienzo et al. 
(2004) 



NGC 


Mrr 




Mmax 




Mo 


Mo 


Mo 


1851 


0.66 


0.666 


0.688 


3201 


0.69 


0.654 


0.706 


4499 


0.70 


0.693 


0.733 


4590 


0.80 


0.704 


0.756 


5272 


0.69 


0.651 


0.698 


5466 


0.74 


0.685 


0.738 


5904 


0.66 


0.627 


0.687 


6362 


0.66 


0.652 


0.686 


6809 


0.69 


0.638 


0.676 


6934 


0.70 


0.673 


0.705 


7078 


0.77 


0.688 


0.767 


7089 


0.66 


0.628 


0.699 



ulation of RR Lyrae variables dominated by ZAHB stars to a 
population dominated by evolved objects. 

Hereafter, we assume that all clusters for which 
f(RR) >0.10 (that is S(RR) >19.7) have an RR Lyrae pop- 
ulation dominated by stars that are still close to their ZAHB 
location. While as mentioned this limit is rather arbitrary, we 
found very similar results with other values in the range 0.04- 
0.10. 

After identifying GCs for which we may assume that most 
RR Lyrae remain close to the ZAHB, we can perform a san- 
ity check of the minimum, median, and maximum masses that 
we adopted for HB stars by comparing them with those derived 
from pulsational properties. Table|9]lists the different mass val- 
ues for the GCs analysed by Di Criscienzo et al. (2004). We 
found that masses derived from pulsational properties are very 
similar to the median masses for the metal-rich clusters (which 
typically have an HB extending on both sides of the instability 
strip), and closer to the maximum masses for metal-poor clus- 
ters, which have predominantly blue HBs. While this result is 
unsurprising because we use the same HB evolutionary models 
as Di Criscienzo et al. (2004), who performed a similar com- 
parison, at least it shows that there are no gross errors in the 
way we applied these models. 

3.2. Discussion 

One of the most intriguing characteristics of HBs is the di- 
chotomy in the distribution of the GCs in terms of the mean 



periods of variables pulsating in the fundamental mode (RRab), 
discovered by Oosterhoff (1944). We note that while the 
Oosterhoff dichotomy applies not only to GCs, but also to field 
stars in our Galaxy (e.g., Szczygiel et al. 2009), populations 
of RR Lyrae sampled in other (smaller) galaxies such as the 
Magellanic Clouds or the dSphs have mean periods that are in- 
termediate between the two Oosterhoff groups (Ooi and Ooii, 
e.g., Pritzl et al. 2004, Catelan 2009). A quite extensive discus- 
sion of various properties of GCs in relation to the Oosterhoff 
dichotomy can be found in Sect. 6 of Catelan (2009). The in- 
formation collectively is quite difficult to interpret, because of 
significant uncertainties in many of the basic parameters con- 
sidered (in particular, ages, but also metallicities), of small 
number statistics, and because average periods within a clus- 
ter depend on a complex combination of factors, mainly re- 
lated to RR Lyrae possibly belonging to different evolution- 
ary phases. In spite of this, many authors have searched for 
a simple explanation of the Oosterhoff dichotomy, and sev- 
eral of them suggested that this might be caused by a com- 
bination of ages and metallicities, the OoII group being at- 
tributable to old, metal-poor GCs, and the Ool to a population 
of young, metal-intermediate clusters (see e.g., Lee & Carney 
1999). We used our data to test whether this interpretation of 
the Oosterhoff dichotomy is consistent with the HB scenario 
we consider in this paper. We first replotted (upper panel of 
Fig. [TST l the metallicity-median mass diagram of Fig. [TO] but 
this time using different symbols for GCs that are either rich 
or poor in RR Lyrae: black and grey circles are for f>0. 10 and 
0.04<f<0.10, respectively. The GCs with f<0.04 are divided 
according to the HBR ratios, (blue) triangles and (red) squares 
representing HBR greater and smaller than 0, respectively. Not 
unexpectedly, RR Lyrae-rich GCs occupy a well defined locus 
in this diagram, roughly represented by a diagonal strip; GCs 
below this strip have blue HBs, while those above the strip have 
red HBs. We note that the anomalous behaviour of NGC 6535, 
being a BHB cluster with very low RR Lyrae fraction and lo- 
cated within the instability strip (shaded box) probably reflects 
a small amount of data in the CMD of this cluster 

As shown by Fig. [TO] this diagram is essentially an age- 
metallicity diagram (with some caveats that are described more 
clearly in the next section). Given the distribution of GCs in age 
and metallicity, there are very few galactic GCs with young 
ages and -2<[Fe/H]<-1.6. This lack of GCs with suitable pa- 
rameters divides the galactic population of RR Lyrae rich GCs 
into two ensembles, which can be identified with the Oosterhoff 
groups, confirming what had already been proposed by many 
other authors. However, ages of both groups have a much larger 
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Fig. 13. Median mass and age as function of metallicity; different symbols in this case are for GCs with different RR Lyrae 
fractions (black: f(RR)>0. 10, grey: 0.04<f(RR)<0. 10) and HBR ratios (blue triangles: HBR>0, red squares: HBR<0). The shaded 
area is the strip containing GCs rich in RR Lyrae. 



spread than given by this simple consideration. This is because 
of two results: (i) the diagonal strip defined by the RR Lyrae- 
rich clusters has a considerable width, much larger than that 
caused simply by the range in colour of the instability strip it- 
self; this indicates that there is an intrinsic spread in the masses 
of ZAHB stars, which is a basic property of HBs that we ex- 
ploit throughout this paper This is important in particular for 
the case of metal-poor clusters, which are able to produce an 
RR Lyrae rich population over a rather wide range of ages, (ii) 
Given the interplay between the original mass and the mass- 
loss variation with metallicity, the strip occupied by RR Lyrae- 
rich clusters overlaps twice with the locus occupied by old 
GCs in the metallicity-median mass diagram. This implies that 
there are old RR Lyrae-rich GCs that are either metal-poor (the 
classical Oo ii clusters), or in a restricted range of metallicity 
around [Fe/H]=-1. Hence, the most metal-rich among the Oo i 
GCs belong to the old population, and are actually as old as 
the most metal-poor Oo ii clusters. This is clearly shown in the 
age-metallicity diagram shown in the bottom panel of Fig. [13] 
where again we plotted with different symbols clusters that are 
either rich or poor in RR Lyrae. 

This description only refers to those GCs that are very rich 
in RR Lyrae. There are many GCs with only a few RR Lyrae: 
as mentioned above, in these cases (generally metal-poor clus- 



ters with predominantly BHB) most variables can be identified 
with stars that have evolved off the ZAHB. In these cases, we 
should expect long mean periods, and the cluster to be identi- 
fied as OoII (e.g., Lee et al. 1990). Furthermore, there are pe- 
culiar cases (including e.g., NGC 6388 and NGC 6441) where 
the RR Lyrae are actually HB stars with masses much lower 
than expected for their metallicity and age (see next section). 
All these cases complicate the simple picture described above, 
as discussed by Catelan (2009). 

To summarise, although the Oosterhoff dichotomy indeed 
depends on the peculiar distribution of galactic GCs in the age- 
metallicity plane, and most of the Oo ii clusters are old while 
most of the Oo i clusters are young, it is not possible to iden- 
tify tout-court the Oosterhoff groups with groups of GCs hav- 
ing different ages (a similar conclusion was drawn by Catelan 
2009). 

4. Ages from median HB colours 

After demonstrating that the median colour of the HB is deter- 
mined mainly by cluster metallicity and age, provided a suit- 
able mass-loss law is adopted, we attempt to derive ages for 
GCs from their median HB colour. This procedure is almost 
identical to that considered by several authors (e.g., Lee et al. 
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Fig. 14. Comparison of ages derived from HB for HST and 
ground-based samples. 

1994), but it is now applied to a large database of GCs in a 
uniform way. 

These ages can be obtained by simply adding to the median 
mass of the HB stars the mass lost during the previous evo- 
lution (Eq. fTTT i. then subtracting the mass difference between 
stars at the tip of the RGB and the TO (Eq. [TOll, and finally in- 
verting the age-metallicity-TO mass relation (Eq.|9]l- Of course, 
since the mean mass-loss law adopted throughout this paper (a 
linear dependence on [Fe/H]) was obtained by comparing me- 
dian masses on the HB with those at the tip of the RGB, where 
mass loss is neglected, on average the two sets of ages should 
agree. However, any difference in the ages derived from the 
HB and RGB for a given GC indicates peculiarities for that 
GC. Taking an extreme view, were not age the main parame- 
ter affecting the HB morphology, we might have no correlation 
(or even an anticorrelation) between MS and HB ages, only 
their mean values being in agreement. The derivation of ages 
from HB is then crucial to understanding the second param- 
eter issue. These ages (expressed in our usual relative scale) 
are listed in Cols. 3 and 4 of Table [TOl for HST and ground- 
based data, respectively. We compared these two age estimates 
in Fig. [14] The two sets of ages agree very well each other, 
with the exception of NGC 2808, for which the ground-based 
HB age is underestimated because the limiting magnitude of 
the relevant photometry is too bright to include the faint end 
of the HB. Once this cluster is dropped, the mean difference 
is 0.003+0.009 (rm.s=0.046); assuming the same error in both 
age estimates, this implies an internal error of 0.032 for each of 
them. In the following, we adopt as our most robust estimates 
of the HB ages those determined from HST data, whenever 
available, otherwise we use values obtained from ground-based 
data. 



0,8 1 1,8 

Age (MS) 

Fig. 15. Comparison of ages from HB and from main sequence 
for both HST and ground based samples; filled squares and cir- 
cles are GCs with My < -8 mag and My > -6.5, respectively, 
while empty triangles are for -8.0 < My <-6.5. Note that the 
scatter around a 1 : 1 correlation is larger for clusters of interme- 
diate luminosity (see text). 

Table 10. Clusters' ages from main sequence fitting and from 
HB both for HST and ground based data-set. Ymed is given 
in Table [12] The complete table is available only in electronic 



form. 


Cluster 


Age(MS) 


Age(HB) 


Age(HB) 


Age(Y) 






HST 


ground 




NGC 104 


0.95 


0.877 


0.856 


0.922 


NGC 288 


0.90 




1.175 


0.958 


NGC 362 


0.80 


0.744 


0.738 


0.789 


IC 1257 




0.880 






NGC 1261 


0.79 


0.735 


0.741 


0.780 



In Fig. [15] we compare these HB ages with those deter- 
mined from the main sequence. The agreement is fairly good: 
the r.m.s. of the difference in ages from HB and from the main 
sequence is 0.094 for HST and 0.105 for ground-based data re- 
spectively: this is approximately equivalent to 1 Gyr While this 
agreement may appear satisfactory, it is actually much larger 
than internal error barsQ 

There are various ways to convince ourselves that this is 
indeed the case. First, we compared in Fig. [16] the HBs of 
two GCs (NGC 1904 and NGC 6934) that have very similar 
metallicities ([Fe/H]=-1.58 and [Fe/H]=-1.56, respectively) 
and relative ages from MS (0.89 and 0.88). It is clear that their 

^ This indicates that the scatter about the mean mass-loss relation 
in Fig.[TT]is indeed real, as mentioned in Sect. l2.4l and not simply due 
to observational errors. 
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Fig. 16. Dereddened HST CMDs for stars on the HB of 
NGC 6934 andNGC 1904. 



HBs are vastly different. The median mass values we obtain are 
0.606±0.003 and 0.671+0.008 for NGC 1904 and NGC 6934, 
respectively. The two clusters also have a very different fre- 
quency of RR Lyrae, i.e., f(RR)=0.02 and 0.38. 

Second, as shown by Fig. [15] the scatter of ages from MS 
and from HB is much larger for GCs of intermediate luminosity 
and smaller for both brighter and fainter clusters. Furthermore, 
there is a clear trend of this scatter with metallicity. This is 
shown in Fig. [T7J where GCs of different metallicity are plot- 
ted separately in the three panels. In this figure, we also plotted 
the typical error bars as well as arrows showing the impact of 
errors in the assumed metallicity and He content. These errors 
take consistently into consideration the impact of the He con- 
tent on the MS colour, the TO point, and the HB. While the 
uncertainties in the assumed metallicity have a relatively mod- 
erate impact on these age estimates, those in the He content 
very significantly affect the age estimates from HB. This ef- 
fect is particularly large at intermediate metallicity, which is 
where we found the largest discrepancies between the two sets 
of ages. This suggests that the He content might be responsible 
for these discrepancies (although we cannot exclude that other 
terms, e.g., the total CNO abundances, might contribute). We 
also note that the errors in the ages from MS and HB due to He 
are anticorrelated with each other: if the true He abundances 
are higher than we assumed, the MS ages are underestimated, 
while the HB ages are overestimated. This result could be ex- 
ploited to simultaneously derive the values of Y^ed and age 
for each cluster, by ensuring agreement between these two age 
derivations. The age values, Age(Y), obtained by applying this 
procedure are listed in the last column of Table[TO] Age(Y) val- 
ues are typically close to MS ages, the largest differences being 



Fig. 17. Comparison of ages from HB and from main sequence 
for the whole sample (HST-Hground based); typical error bars 
are shown at bottom right corner of each panel. The thick 
(black) arrows show the impact of a 0.02 variation in the as- 
sumed He content, while the thin (red) ones represent the con- 
sequence of a 0.1 dex error in [Fe/H]. 
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Fig. 18. Age-metallicity relation for different groups of glob- 
ular clusters: outer halo clusters (green pentagons), inner 
halo clusters (blue circles), disk/bulge clusters (red squares). 
Magenta circles are GCs associated with dSphs. Different sym- 
bol sizes are used for clusters of different luminosity. 
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Table 11. R' parameter as derived from both HST and ground 
based data (note the different definition with respect to the 
usual R parameter). The final adopted value <R'> is the 
weighted average between the two estimates, when available. 
Y(<R'>) is the corresponding He abundance. The complete ta- 
ble is available only in electronic form. 



NGC /IC 


R'(HST) 


R' (ground based) 


<R'> 


Y(<R'>) 


104 


0.69±0.07 


0.68+0.09 


0.69+0.07 


0.253+0.012 


288 




0.69+0.14 


0.69+0.14 


0.250+0.024 


362 


0.75+0.08 


0.63+0.15 


0.75+0.08 


0.260+0.012 


1257 


0.49+0.14 




0.49+0.14 


0.213+0.024 


1261 


0.62+0.10 


0.69+0.09 


0.62+0.10 


0.238+0.017 



for a few clusters with very blue HBs, such as NGC 6254 and 
NGC 1904. 

In principle, errors in Age(Y) should be obtained by 
quadratically summing the internal errors in MS and HB ages; 
the mean quadratic value determined in this way is 0.063. To 
confirm the reliability of this uncertainty estimate, we plot the 
age-metallicity relation for different groups of GCs in Fig [18] 
outer halo GCs, inner halo GCs, disk/bulge GCs, and GCs as- 
sociated with dSphs (for the classification of the GCs in these 
different populations see Carretta et al. 2009d). Briefly, we 
replot the same age-metallicity diagram, as previously done 
in Carretta et al. (2009d) but with the new age values, i.e. 
Age(Y): we found that for the disk/bulge GCs the relation 
becomes tighter, with an r.m.s around a least squares fit of 
0.036 with respect to 0.042 previously derived (note that as- 
suming the age values directly retrieved from Marm-Franch et 
al. 2009, the scatter around the fit line is ~0.050). This im- 
provement also supports the identification of this third parame- 
ter with variations in the median He abundances. Furthermore, 
the very small dispersion, well within the above cited error of 
0.063, suggests that the observational uncertainties in the ages 
from MS's were most likely overestimated. We consider these 
Age(Y) values to be the most robust estimates of the relative 
ages currently available: they incorporate updated photometric 
data and metallicities, and take into account cluster-to-cluster 
variations in the He content. We however postpone to a future 
paper a full discussion of the age-metallicity relation for the 
different Galactic populations of GCs. 

5. He from HB and R-parameter 

Are there other arguments supporting the idea that the third 
parameter required to explain the median colours and masses 
of HB stars is indeed the He content? Is there any available 
evidence on variations in the median He content between dif- 
ferent GCs? The most classical derivation of He abundances 
for GCs is obtained using the so-called R-parameter method 
(Iben 1968). The R-parameter is the ratio of the number of HB 
stars, Nhb, to the number of RGB stars brighter than the HB in 
MboU Nrgb- This ratio depends on the He content, because the 
higher the He content, the brighter the HB, the smaller the num- 
ber of RGB stars brighter than this level, and hence the larger 
the value of R. The R-parameter can be calibrated using mod- 
els that include the variation in the lifetime of RGB and HB 



stars with metallicity. The most updated calibration of R is that 
by Cassisi et al. (2003). This calibration was used by Salaris et 
al. (2004) to produce R values and average He abundances for 
57 GCs based on the same HST photometric database consid- 
ered here. Salaris et al. (2004) set an upper limit to the r.m.s. 
scatter of Y among these GCs of 0.019. By itself, this upper 
limit to the rm.s. scatter is fully consistent with the rm.s. scat- 
ter of 0.014 for Fmed that we derived by matching HB and MS 
ages. It would be rewarding to identify some degree of correla- 
tion between the Y values of individual GCs derived by the two 
techniques. However, the correlation is actually poor, probably 
because the Y values, derived using the R-method, have large 
error bars. This is mainly due to limited statistics, though in a 
few cases other error sources (contamination by field stars, dif- 
ferential reddening, photometric limits) contribute to the noise 
in Y determinations. 

We note that the statistical error in the R parameter can in 
principle be reduced by adopting a different definition. With the 
usual definition, on average R~ 1.5, there are more stars on the 
HB than on the RGB brighter than the HB, which is then a ma- 
jor contributor to the random noise. However, we may use a dif- 
ferent definition of R, where the error caused by the small num- 
ber of stars on the RGB is substantially reduced, but has a sim- 
ilar dependence on the helium content than the R-parameter. In 
practice, we may define R' = Nhb/N^q^, where N^q^ is the 
number of stars on the RGB brighter than V(HB)h-1. Roughly, 
N^QQ is twice the value of A^rgb, so that R' - Nhb /-^rgb ~ 
It is then clear that the statistical errors in R' are smaller than 
in R. Furthermore, checks of the stellar models of Bertelli et 
al. (2008) indicate that R' has roughly the same (fractional) de- 
pendence on Y as R, with only one important modification: the 
RGB bump is brighter than V(HB)-i-l in all GCs considered in 
this paper Hence, R' has a simpler (roughly linear) dependence 
on [Fe/H], unlike R. We then decided to repeat the test using Y 
values inferred from R', rather than R. 

Table [TT] gives the values of R' deduced from HST and 
ground-based data; since completeness is more of a concern 
for this second set of photometry, we give preference to the 
HST values, whenever possible. When making this derivation, 
we simply used the HB level listed by Harris (1996). The er- 
rors related to individual entries of Table [TT] are simply those 
from statistics, and do not include the effects of uncertainties 
in decontamination by field stars, differential reddening, and 
photometric limits. In a few cases, these effects are known and 
can be taken into consideration. On this basis, we exclude from 
our estimates the ground-based observations for NGC 2808, 
NGC 5986, and NGC 6266, and the HST data for NGC 2419, 
where the photometry is not deep enough to reach the ex- 
treme hot/faint tail of the HB. For NGC 5927, the ground-based 
photometry is also not good enough to discriminate between 
HB and RGB stars. The R' values are underestimated in these 
cases, and we did not consider them in our discussion. 

In Fig. [19] we compare the average values of R' with the 
values of R obtained by Salaris et al. (2004). There is a good 
correlation between these two sets of data, which improves if 
we consider the offset between metal-rich and metal-poor clus- 
ters caused by the luminosity of the RGB bump; however, we 
recall that this new determination is not completely indepen- 
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Fig. 19. The R' values plotted against the R values from Salaris 
et al. (2004); filled and empty squares are for GCs with [Fe/H] 
larger and smaller than -1.0 dex, respectively. Note the differ- 
ent definitions of R' and R (see text), 



Fig. 20. The R' values plotted against [Fe/H]; filled and open 
symbols are GCs with or without age estimate, respectively. 
The solid lines correspond to Y=0.225, 0.250, 0.275. 



dent of the Salaris et al. one, since both use the HST snap- 
shot photometry. A reasonable good correlation is also obtained 
with the older determinations of R by Sandquist et al. (2000) 
and Zoccali et al. (2000). 

The last column of Table [TT| gives the value of Y that we 
could derive from these R'. These were obtained by scaling 
the calibration of Cassisi et al. (2003) for the ratio of N^q^ to 
A^RGB, and taking into account that the RGB bump is brighter 
than V(HB)h-1 in all GCs, while being fainter than V(HB) for 
GCs with [Fe/H]> -1.1. In Fig.|20]we have plotted the values 
of R' against the metallicity [Fe/H], and compared them with 
the calibrations we used. We note the absence of an obvious 

bump at [Fe/H] 1.1, which is present in the case of the R- 

parameter. 

In Fig. [21] we plotted the values of Y{R') obtained follow- 
ing this procedure against the median He abundance Fmed re- 
quired by placing ages estimated from the HB in agreement 
with those determined from the MS. There is some correlation 
between the two determinations of Y (linear correlation coef- 
ficient r=0.23 over 64 GCs, which is significant at a level of 
confidence of > 95%). The scatter is clearly not negligible, as 
expected based on the rather large errors in Y(R'). However, 
we note that the meaning of these two quantities is not neces- 
sarily the same. To show this, we note that the strongest con- 
tribution to a test performed by comparing the two sets of 
He determinations is from NGC 6388 and NGC 6441; in both 
these clusters, Y(R') is far larger than Ymed than the (statisti- 
cal) errors. The large value of Y(R') is explained by HB at the 
RR Lyrae colour being very bright in these clusters (and indeed 
the RR Lyrae have anomalously long periods: Pritzl et al. 2002, 




0,15 I ' L_ 
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Fig. 21. The values of Ymed are plotted against Y(R'); filled 
and open squares are for errors smaller and larger than 0.02, 
respectively. 



2003). On the other hand, the RR Lyrae are most likely in the 
tail of the distribution of high He abundances in these clusters, 
whose HB is predominantly red. In these cases, Y(R') is then 
expected to be larger than Fmed, as is indeed observed. We also 
note that a much stronger correlation exists between Y(R') and 
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Fig. 22. Observed CMDs for NGC 4833 and NGC 5024 (upper 
panels) and synthetic diagrams lower panels), obtained with a 
Gaussian-like distribution mass loss with standard deviation of 
0.02 Mq and 0.07 Mq, respectively. 

the average of Ymed and Ymax- In this case, the correlation co- 
efficient is r=0.294 for the same 64 GCs, which is significant at 
a level of confidence of ~ 99%. 

We conclude that the hypothesis that the He abundances ob- 
tained by matching ages from MS and HB is compatible with 
current estimates of the He abundances from the R-method, but 
that more extensive data sets and more careful comparisons 
(taking into consideration other possible indication of the He 
abundance, see e.g., Catelan et al. 2009) are needed to defini- 
tively settle this issue. 

6. The colour and mass spread along the HB 

6.1. Observational data 

The previous discussion suggests that at least three parameters 
(metallicity, age, and most likely helium content) are required 
to provide a detailed explanation of the median colours of the 
HBs. We intend to demonstrate that this selection of parameters 
successfully explains the extension of the HBs. 

To illustrate our approach, in Fig. |22] we compare the BV 
colour magnitude diagrams of NGC 4833 and NGC 5024. 
These two clusters have very similar ages (1.01 and 1.04, re- 
spectively) and metallicity ([Fe/H]=-1.89 and -2.06 respec- 
tively), similar values of both median colour and mass of stars 
on the HB (0.655 versus 0.658 Mq), but yet very different HBs. 
The main difference between these two clusters is in the min- 
imum mass of stars along the HB (0.561 vs. 0.649 Mq). We 
can explain the short HB of NGC 5024 in terms of a small 
star-to-star scatter in mass loss (see e.g.. Rood et al. 1973, and 
many others since). The bottom left panel of Fig. |22] shows a 




B-V 



Fig. 23. Cumulative distribution of (B-V) colours from 
colour-magnitude (dashed) and synthetic (solid) diagrams for 
NGC 5024 and NGC 4833; the resulting probabiHties from 
Kolmogorov-Smirnov (KS) test are also shown. 



synthetic HB obtained by assuming that the mass lost during 
the RGB phase is distributed like a Gaussian, with a standard 
deviation of 0.02 M© (this is the minimum value required to 
reproduce the shortest HBs). While this synthetic HB matches 
well the one observed in NGC 5024, it is clearly a poor match 
to that of NGC 4833, because it fails to describe its very ex- 
tended faint tail. The difference between the two clusters is 
not simply in terms of the spread of mass loss. This is shown 
in the bottom right panel of Fig. [22] where we plotted a syn- 
thetic HB obtained assuming a mass loss spread of 0.07 Mq, 
required for an HB extending as faint as the observed one for 
NGC 4833. However, this synthetic HB has by far too many red 
stars, as shown by a comparison of the cumulative distributions 
in colours (see Fig.l23]lri A Kolmogorov-Smirnov test shows 
that the observed colour distribution disagrees with that of the 
synthetic HB (at about 98% level of confidence). In contrast, a 
much better agreement is obtained between the observed and 
synthetic HBs for NGC 5024. We may conclude that at least in 
the case of the pair NGC 4833-NGC 5024, a third parameter 
other than age and metallicity is required, and that this param- 
eter mainly affects the minimum mass along the HB. 



* Small excesses of red stars along the HB - much smaller than 
found for NGC 4833 - are found often in our comparisons with syn- 
thetic HBs. This same problem was found by other similar analysis 
(see the discussion in Castellani et al. 2005). It might indicate either a 
deficiency of the adopted HB models, or a non-Gaussian distribution 
of the random terms of the mass loss, perhaps because it would be 
more appropriate to consider a Gaussian distribution in the mass-loss 
efficiency parameter (D'Cruz et al. 1996). 
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Fig. 24. Same as Fig. |22l but for the GCs NGC 6171 and Fig. 26. Absolute magnitudes (My) plotted against the difFer- 
NGC 6723. The synthetic diagrams (lower panels) were ob- ence between the median and the minimum mass along the HB 
tained with a Gaussian-like distribution mass loss with standard (<5M). In this case, we also separate GCs older (red triangles) 



deviation of 0.03 and 0.07 Mq, respectively. 



and younger (blue filled dots) than 0.92. 




again vastly diff'erent HB, due to a very different value for the 
minimum mass (0.629 vs. 0.558 Mq). We can reproduce the 
observed HB of NGC 6171 by a synthetic HB with a small 
Gaussian spread of 0.03 in the mass lost (see Fig.l24li.The 
shape of the colour distribution of NGC 6171 is not perfectly 
reproduced, again because the synthetic diagram predicts too 
many red stars, as shown by Fig.|25j however, given the small 
size of the sample, the result of the Kolmogorov-Smirnov test 
returns a 32.7% probability that observed and synthetic HBs 
are drawn from the same population, which is not significant. 
On the other hand, the distribution of colours of HB stars along 
the HB of NGC 6723 is clearly bimodal (see also Rood & 
Crocker 1985), and there is no way of reproducing it assum- 
ing a Gaussian distribution of mass lost. In this case, the re- 
sult of the Kolmogorov-Smirnov test is that the probability 
that a Gaussian distribution (with a dispersion of 0.07 Mq) 
reproduces observations is 0.2%. Hence, the case of the pair 



NGC 4833-NGC 5024 is not alone. 



Fig. 25. As for Fig.|23] cumulative distributions and KS values 
for NGC 6171 and NGC 6723. 



A similar comparison can be made using more metal-rich 
clusters. For this purpose, we selected the pair NGC 6171- 
NGC 6723, which also have very similar metallicities and ages 
([Fe/H]=-1.03 and -1.10, Age=0.99 and 1.01, respectively), 
similar median mass on the HB (0.650 and 0.644 Mq), but 



6.2. The correlation between mass spread and 
absolute magnitude 

We propose that the cases discussed in the previous subsection, 
extracted from several decades of discussions, convincingly 
show the need for a third parameter in addition to metallicity 
and age to explain the blue extension of the HBs. Hereinafter, 
we quantify the spread in mass loss by the difference between 
the minimum and median mass of stars along the HB, and call 
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this quantity 6M 0. We note that 6M has been estimated by 
assuming a constant He abundance, hence it cannot be used 
directly to derive the spread in He within the cluster 

It has long been known that GCs with a large spread of 
masses along the HB are bright and massive (Fusi Pecci et al. 
1993). Recio-Blanco et al. (2006) obtained a good correlation 
between the highest effective temperature along the HB and 
the cluster absolute magnitude My. Given that the highest ef- 
fective temperature along the HB is well (anti-)correlated with 
the minimum mass of HB stars, we should expect a good anti- 
correlation between 6M and My. In Fig.|26l we show the run 
of (5M with My (My values were taken directly from the Harris 
1996 catalogue), for both the HST and ground-based data sets. 
There is indeed a very significant anti-correlation between 
these two quantities: large values of 6M are only found among 
bright clusters. The correlation is much cleaner if we restrict 
ourselves only to those GCs with age estimates, which have 
the highest quality colour-magnitude diagrams. Hereinafter, we 
use these GCs only, and combined ground-based and HST data. 
However, by comparing the two data sets, we noted that the 
anti-correlation between 6M and My is stronger when HST 
data are used (r=0.68 for 45 clusters) than if ground-based data 
are considered (r=0.54 for 46 clusters). This stronger corre- 
lation occurs because the extreme BHB stars are sometimes 
below the detection limit of the ground-based observations. 
Hence, when considering the two samples (whenever possible), 
we give preference to the HST data. 

For this whole sample of 65 GCs (making up almost half of 
the total number of Galactic GCs), the linear relation between 
6M and My is 

5M = -(0.102 ±0.026) -(0.020 + 0.003) My (12) 

with an extremely significant correlation coefficient of r=0.63. 
The r.m.s. scatter of individual points along this mean relation 
(0.026 Mq) is not much larger than expected from internal er- 
rors, suggesting that it is mainly due to observational errors. 
However, part of the scatter is certainly real. We note that there 
are a few clusters (namely NGC 104=47 Tuc, NGC 5024; and 
to a lesser extent NGC 362 and NGC 5272=M 3) that have a 
small spread of masses along the HB in spite of there being 
quite massive; one of them has indeed been used in the com- 
parisons of the previous subsection. There is little doubt that 
the spread in mass along the HB is small, at least for the two 
most extreme cases. We conclude that the additional parameter 
required to explain the extension of the HB, while quite closely 
related to the overall cluster luminosity, is actually a separate 
one. 

6.3. He variations required explaining tlie HB widtli 

In the remainder of this section, we examine the possibility that 
the additional parameter (to both metalUcity and age) determin- 
ing the spread in mass for stars along the HB is the variation in 

' In principle, Mn,in should give a more correct estimate of the 
spread in mass along the HB. However, we preferred to use M,^e(] in 
our discussion for the reasons explained in Sect. 2.3, even though in 
some cases it may systematically underestimate the mass correspond- 
ing to stars with primordial He abundances 



the He abundance, related to various stellar generations in GCs, 
which combines with a small (~ 0.02 Mq) random spread in 
mass loss. This last value was adopted because it is roughly the 
value required to explain those GCs with the minimum spread 
in mass for stars along the HB. In principle, this value may vary 
from cluster to cluster, depending e.g., on cluster concentration 
(see e.g., Fusi Pecci et al. 1993). However, we wish at present 
to keep our assumptions to a minimum. 

To test the hypothesis presented above, we first derived the 
spread in He required to explain the observed spread in colours 
and masses along the HB. In the next section, we discuss the 
evidence provided by chemical abundances that might support 
this hypothesis. When deriving the spread in He, we should 
take into account that variations in He abundances have impor- 
tant effects on our analysis. 

The first step of our procedure is quite simple. Since we as- 
sumed that an intrinsic spread in the mass loss equal to 0.02 Mq 
is a universal phenomenon, we corrected the observed mass 
spread for this effect by subtracting this value in quadrature 
from the observed 6M. In those few cases where 6M< 0.02, 
we simply assumed that the corrected value is 0. The coiTected 
spread is then attributed to a variation in the He abundance, 
assuming that the same mass-loss law is applicable to all stars. 

However, to derive the variation in He, we should take into 
account that the masses derived for HB stars are themselves 
functions of the adopted He abundance. Unfortunately, the Pisa 
evolutionary tracks were computed for a single (not constant) 
value of the He abundance for each value of the metal abun- 
dance, and cannot then be used to estimate this correction. We 
therefore used instead the isochrones by Bertelli et al. (2008), 
which while assuming a unique value for the mass loss along 
the RGB (preventing its use to derive the relation between 
masses and colours along the HB), do however provide data 
for different He abundances. We then combined the two sets of 
models to produce the following correction formula for AM 

AM(Y) = AM-AY(1.976-Hl.982[Fe/H]-H0.562[Fe/Hft.(13) 

In this equation, AY is the variation in He abundance with re- 
spect to the reference value (Y=0.25, roughly the cosmolog- 
ical value). In this case, we assume that all GCs started their 
evolution with an He abundance close to that resulting from 
the big bang nucleosynthesis (Y^O.25, WMAP, Spergel et al. 
2007 -see their Table 4), but that there might be significant 
star-to-star variations even within the same GC. The procedure 
for estimating this He abundance is described in Sect. |5] Of 
course, the primordial He abundance in some clusters may ac- 
tually be higher than the cosmological one. This is supported 
by both simple arguments based on a rather constant AF/AZ 
value throughout the Galactic evolution (Chiosi & Matteucci 
1982; Balser 2006; Casagrande et al. 2007), and indirect mea- 
surements in GC stars, e.g., those obtained using the R-method 
(see Sect. |5]l. However, these variations in primordial He are 
expected to be quite small (a few hundredths in Y), so that this 
assumption is not critical. In turn, AM(Y) may be derived from 
the evolutionary mass at the tip of the RGB, corrected for the 
uniform mass-loss law along the RGB. In practice, we use the 
equation 

AY = AM(Y)/(1.245 + 30.747 ZqIO^^""™), (14) 
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Fig. 27. Absolute magnitudes vs. Ymax for HST and ground- 
based samples; as in all previous figures, clusters are separated 
according to their ages. 



Table 12. Spread in He and mass required explaining the HB. 
The complete table is available only in electronic form. 



Cluster 


SM(Y) 


-'- med 


Y 

^ max 


dY 


6(B-V) 


6{V-I) 




Mo 








mag 


mag 


NGC 104 


0.000 


0.234 


0.234 


0.000 


0.000 


0.000 


NGC 288 


0.016 


0.280 


0.292 


0.012 


0.005 


0.007 


NGC 362 


0.059 


0.243 


0.289 


0.046 


0.018 


0.025 


NGC 1261 


0.068 


0.244 


0.297 


0.053 


0.021 


0.029 


NGC 1851 


0.063 


0.247 


0.295 


0.048 


0.019 


0.027 



where Z© is the solar metal abundance, and [m/H] is the over- 
all metallicity, which we assumed to be [m/H]=[Fe/H]-i-0.28, 
the second additive term taking into account the overall en- 
hancement of Qf-elements in GCs (see Carretta et al. 2009d). 
These two equations can be solved recursively, convergence be- 
ing quite fast. 

Table [12] lists the values of Y^ed, Y^ax, and SY-YmiLx- 
Ymed obtained by this procedure for each cluster. We recall 
that in some cases Y„ied values might overestimate the mini- 
mum Y in the cluster, and then underestimate SY. The highest 
He abundances found using this procedure is Y„iax=0.347 for 
NGC 6273. The value we find for NGC 2808 is Y^^„^0333; 
this value is discussed in the following section. 

In Fig.|27] we plotted the Ymax values against the absolute 
magnitude My. As expected, the correlation is good, yielding 

Y^ax = (0.180 + 0.025) - (0.0134 + 0.030)My (15) 

with a very significant correlation coefficient of r=0.49 for 65 
clusters. 



6.4. Colour spread along the main sequence 

The spread in helium required to explain the spread in masses 
along the HB should cause a broadening of the MS. We esti- 
mated this expected colour spread in the MS at absolute mag- 
nitude My - +8, by using the following formulae we derived 
by fitting data from the Bertelli et al. (2008) isochrones 

6(B - y) = 57(0.350 - 0.035[Fe/H]) (16) 

and 

6(V - /) = 57(0.672 + 0.099[Fe/H]). (17) 

The colour spreads in Z? - V and V -I for each cluster are given 
in Cols. 6 and 7 of Table [T2l respectively. Again, we recall that 
in some cases these spreads may be underestimated, because 
6Y may itself be underestimated. This spread is very small, be- 
low detectability, for most of the GCs. The largest spread is 
expected for NGC 6273 (0.035 mag inB-V and 0.043 mag in 
V-I); this cluster is affected by a strong differential reddening, 
which complicates the detection of this spread. 

Most GCs have not yet been scrutinised in enough detail, 
but in a few cases we may compare these predictions with 
observations. For instance, in the case of NGC 2808, Piotto 
et al. (2007) found a spread of ~ 0.1 mag in the F475W - 
F814W colour from very high quality ACS data. Since we ex- 
pect that 6(F475W-F8 14W)/6(V-I) ~ 1 .5, the spread we pre- 
dict from our analysis of the HB is d(F475W-FSUW) ~ 0.05, 
which is roughly half the spread indeed observed. Part of this 
difference can be attributed to the median colour of the HB of 
NGC 2808 (B-Vmed = 0.024) not corresponding to the red HB, 
but rather to the BHB, and in our framework is then interpreted 
as a moderately He enriched population (Y=0.273). The He- 
poor population is however present, corresponding to the red 
HB, which makes up almost 40% of the cluster HB population. 
The total spread in colour along the MS is then expected to be 
larger than given simply by 6Y = F^ax - ^med; a more appropri- 
ate estimate of the expected colour spread is 6(V - I) ~ 0.046 
and 6(F475W - F814W) ~ 0.07. While this is still somewhat 
smaller than observed, the discrepancy is now small, and could 
be justified by some additional source of scatter (differential 
reddening, binaries, photometric errors) for MS colours or an 
incorrect calibration of the HB. 

Anderson et al. (2009) found that some spread in colour 
is also present in 47 Tuc, from a very comprehensive analysis 
of extensive archive HST data. This spread is small, roughly 
6(F616W - F814W) ~ 0.02, and could only be identified 
thanks to the exceptional quality of the data, and the very care- 
ful procedures used to analyse them. However, the HB of 47 
Tuc is very shorQ so that our estimate for the He spread is 
only 6Y - 0.0 (based on the assumption that the small resid- 
ual spread in masses could be explained by random star-to-star 

47 Tuc actually contains a few very blue HB stars (Moehler et 
al. 2000); they make up ~ 1% of the HB stars in the HST colour- 
magnitude diagram we are considering. These stars are so rare that 
they do not affect our definition of M^;„, which excludes the lowest 
5% of the distribution. Indeed, this case underlines that there most 
likely are stars at the extremes of the HB whose origin is not related 
to their extreme values of the He abundance. 
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variations in the mass loss adopted throughout this paper), and 
we then expect no widening of the MS due to He. A careful 
study of the HB of 47 Tuc by Di Criscienzo et al. (2009) is work 
in progress; early results suggest that it may be more accurately 
represented by assuming that there is a very small spread in Y 
(dY « 0.02). By itself, this is not enough to justify the spread 
in colours of the MS. However, Di Criscienzo et al. also found 
that variations in the total CNO abundances might possibly ex- 
plain both the small spread in colour of the MS and the far 
more evident split in the SGB, also found by Anderson et al. 
(2009; see however Bergbush & Stetson 2009 for a different 
interpretation of this observation). While this comparison on 
the whole supports our result of a small spread in He in 47 Tuc, 
it also indicates that more careful examinations using synthe- 
sis of populations, taking into account the distribution in both 
colours and magnitudes of the stars on the HB, may provide 
additional important information about the properties of GCs. 
Extension of this careful analysis to many other GCs would be 
very helpful. 

7. He and GC chemistry 

7.1. Light elements anti-correlations 

Additional evidence that the spread in colours along the HB of 
GCs is caused by variations in the He abundances can be ob- 
tained by considering correlations with similar spreads in the 
abundances of light elements such as Na and O, Mg and Al. As 
mentioned in the introduction, these spreads are caused by dif- 
ferent generations of stars in GCs, the ejecta from the earliest 
one having polluted the material from which the later genera- 
tion(s) of stars formed. While p-captures at high temperatures 
are clearly required to produce Na and Al, and destroy O and 
Mg, we do not yet have a satisfactory detailed model, and even 
its astrophysical basis is currently debated (either fast rotating 
massive stars or massive AGB stars, experiencing hot bottom 
burning, or massive binaries: see Decressin et al. 2007; Ventura 
et al. 2001; De Mink et al. 2009). We are unable to derive the 
exact mass range of the polluting stars. This is an important 
concern, because the production of He has probably a different 
dependence on the mass of the polluters than both the produc- 
tion of Na and the destruction of O. For instance, if we consider 
the massive AGB scenario. He is mainly produced in the pre- 
vious MS phase, and it is far more abundant in the ejecta of 
the most massive polluters (mass > 5 Mq). On the other hand, 
a very significant production of Na and depletion of O can be 
obtained even within stars of lower mass. Hence, the ratio of 
He to Na production might change from cluster-to-cluster, pro- 
vided that the mass range of polluters changes. With this caveat 
in mind, we can then examine current evidence. 

In principle, the original He content should be attainable 
from direct measurements for each HB star However, this da- 
tum is neither available nor can be easily obtained, save possi- 
bly for a restricted temperature range at ~ 10000 K (Villanova 
et al. 2009; He is heavily depleted by sedimentation in warmer 
HB stars: see e.g., Behr et al. 2000b, Behr 2003). Furthermore, 
even O and Na abundances are not available for most HB stars 
(for a possible exception, see again Villanova et al. 2009), al- 
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Fig. 28. SY as function of IQR for [Na/O] and [Al/Mg], respec- 
tively. Lower limits are also indicated. 

though in this last case abundances could in principle be ob- 
tained for stars on the RHB, within the instability strip, and 
on the coolest part of the BHB. Because of this shortage of 
data, we use values provided by RGB stars. Of course, the im- 
plication is that only statistical properties of the distributions, 
not individual values can be studied. In practice, we only con- 
sider extreme values, with the assumption, consistent with our 
approach, that the bluest (i.e., less massive) HB stars are the 
progeny of the most He-rich TO-stars, and the reddest ones of 
the most He-poor ones; in a future paper, we plan to consider 
in more detail other characteristics of these distributions. 

In the present discussion, we considered a subset of the 
GCs with extensive enough data on the Na-O and Mg-Al an- 
ticorrelations. Most of the GCs considered here are from the 
very extensive study by Carretta et al. (2009a, 2009b); we 
complement this data set with a few GCs from the literature 
(Shetrone & Keane 2000 for NGC 362; Sneden et al. 2004 
and Cohen & Melendez 2005 for NGC 5272 and NGC 6205; 
Sneden et al. 1 99 1, 2000 for NGC 6341; Marino et al. 2008 for 
NGC 6121; Marino et al. 2009 for NGC 6656; and Yong et al. 
2005 for NGC 6752). The relevant data are given in Table □ 
We note that we prefer to use the inter-quartile IQR of the 
distributions (either directly taken from the literature, or ob- 
tained from abundances for individual star^3)^ rather than the 
corresponding values for 90% of the distribution, as we did 
for the extension of the HB. The IQR is indeed a more ro- 
bust indicator, which is less sensitive to small number statis- 
tics and to many upper limits to the abundance determinations. 
This last problem remains an important concern. In general, 
we may consider the IQR determinations as lower limits for 

" We thank A.F. Marino for providing the unpublished data for the 
individual stars in NGC 6656. 
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Fig. 29. Run of [Na/Fe]n,i„ with [Fe/H] for the GCs of our sam- 
ple. The thick superimposed line represents an average line fit- 
ting field stars (see text) 

those GCs with metallicity [Fe/H]< -1.5. Figure |28] illustrates 
the correlations between the spread in He abundances obtained 
from the colour spread of the HB, and the IQR values for both 
the [Na/O] and the [Al/Mg] anticorrelations. As can be seen, 
fairly tight correlations exist, the strongest being between 6Y 
and IQR([Na/0]), which are based on more extensive data sets 
than the IQR([A1/Mg]). This correlation strongly supports the 
current interpretation that the extent of the HB is determined 
by the spread in He content within each GC. However, we 
note that at least for 5Y-IQR([Na/0], the relation has an offset, 
where 6Y is significantly larger than only for clusters with 
IQR([Na/0]> 0.6). 

7.2. Comparison with a dilution model 

We may compare these results with the prediction of a simple 
universal mechanism for the production of the Na-O correla- 
tion, that has only one class of polluters, and for which the 
composition of second generation stars differ only in terms of 
the different dilution of polluted gas by primordial material. 
A similar dilution model has been successfully used to ex- 
plain many features of the Na-O anti-correlation (see discus- 
sion in Prantzos & Charbonnel 2006), including its shape and 
the residual observed Li abundances in stars severely depleted 
in O in NGC 6752 (Pasquini et al. 2005). Once the composi- 
tions of the pristine and processed material are set (e.g., by the 
extremes of the observed distributions), the dilution factor may 
be determined for each star (either from O or Na abundances), 
and the helium production can be inferred, save for a constant 
factor 

In this model, the logarithmic abundance of an element [X] 
for a given dilution factor dil is given by 

[X] = log [(1 - dil) 10'^°^ + dil lO'^pl], (18) 

where [Xo] and [Xp] are the logarithmic abundances of the 
element in the original and processed material, respectively. 
Practically speaking, we may assume (see Carretta et al. 2009a, 
Fig. 17) that the original composition is [Na/Fe]=-0.28 if 
[Fe/H]< -1.5, [Na/Fe]=0 if [Fe/H]> -0.8, and [Na/Fe]=- 
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Fig. 30. Run of [Na/Fe] vs. [O/Fe] expected for different values 
of dil from our dilution model. The insert shows the run of 
[Na/Fe] with [O/Fe] for stars in our FLAMES survey (Carretta 
et al. 2009a). 



0.28+0.4([Fe/H]-i-1.5) if -1.5<[Fe/H]<-0.8; and [O/Fe]=0.5. 
For the processed material we assumed [Na/Fe]=0.6 and 
[0/Fe]=-1.5. We note that while for O we may assume a 
uniform [O/Fe] value, because cluster-to-cluster variations of 
maximum values are small among the GCs of our sample, we 
adopted an original [Na/Fe] ratio that is a function of metal 
abundance, to reproduce both our data and the Na abundances 
observed among field stars (see Fig.l29ll. 

In Fig. [30l we compare the expected run of [O/Fe] vs. 
[Na/Fe] according to this dilution model, with observations for 
the ensemble of stars observed within our project (Carretta et 
al. 2009a). It is important to notice that within this approach the 
production of Na (and of He) saturates at large values of the di- 
lution parameter dil. This has two important implications: (i) 
most reliable estimate of the maximum value of dil is given 
by the O abundances; and (ii) that similar He abundances are 
obtained for a rather wide range of [O/Fe] . On the other hand, 
the minimum value of dil is most accurately given by Na abun- 
dances, because similar values of [O/Fe] are obtained for small 
values of dil. Minimum and maximum values of O and Na 
abundances measured in each cluster are given in Cols. 2, 3 and 
4 and 5 of Table [13] We note that while Na and O maximum 
abundances listed in this Table are somewhat different from the 
values adopted in our dilution model, this has no consequences 
on our discussion because minimum and maximum values of 
the dil parameter (Cols. 6 and 7 of Table [T3]) are obtained from 
the minimum abundances of Na and O, respectively. Both min- 
imum and maximum values of dil contain some uncertainty, 
because only upper limits to O and Na abundances can be ob- 
tained in extreme cases; this caveat is stronger for metal-poor 
stars, and for the smallest GCs, where most of the observed 
stars are quite warm. For instance, we suspect that the mini- 
mum Na abundance is overestimated in the case of M 15, which 
is a very metal-poor GC, and that for this reason we underes- 
timate the 6 dil range appropriate for this cluster Finally, the 
differences between the maximum and the minimum dilution 
are given in the last column of Table [13] Within this model, 
these values should be proportional to 6Y. 
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Table 13. Minimum and maximum values for [Na/Fe] and 
[O/Fe] (see text for references). The dilution fraction is also 
reported. 
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We note that while this simple dilution model predicts quite 
uniform values of 5 dil, and then 6Y, the values we derived 
from the HB exhibit large variations from cluster-to-cluster As 
a consequence, the correlation for individual clusters is poor 
We may also remind that this simple model, which uses a uni- 
versal polluter, also predicts that the Na-O and Mg-Al anticor- 
relations should closely resemble each other, which is not ob- 
served (see Carretta et al. 2009b). 

We conclude that while a second generation of stars pol- 
luted by some stars of an earlier generation exists in all GCs, 
the composition of the polluters is not universal. In some clus- 
ters (such as NGC 2808), the polluters produce large amounts 
of He (most likely with similarly large amounts of Al), and very 
efficiently destroy O and IVlg; while in others (like NGC 6121) 
practically no fresh He (and Al) is present in the material that 
produced the second generation stars. On the other hand, all 
these polluters produce similar amounts of Na. Independent of 
the polluter, we should be able to reproduce these observations, 
within a viable scenario for cluster formation and early evolu- 
tion. 



7.3. The impact of cluster luminosity 

The previous discussion suggests that the spread in the He con- 
tent within stars in GCs are the most important factor determin- 
ing the extension of the HBs. It is important however to under- 
stand the basic physical property of GCs causing this spread 
in He. The existence of a good correlation between the maxi- 
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mum temperature of the HB and the cluster luminosity (Recio- 
Blanco et al. 2006) suggests that the total cluster mass deter- 
mines its ability to have stars of very different He contents, 
assuming that the present cluster luminosity is a good proxy 
for its original mass. In Fig.[3T|we show the correlation that 
we find between the spread in He abundance (as represented 
by 6Y) and cluster luminosity My, here a proxy for the total 
mass. There is a clear strong trend of an increase in 6Y with 
decreasing My (that is an increase in luminosity). The linear 
correlation coefficient is 0.56 for a total of 65 clusters, which is 
significant at a level of more than 99.99 %. The mean regres- 
sion line is 



6Y = -(0.0127 + 0.0024)/Wv - (0.069 + 0.020). 



(19) 



Although the correlation between the spread in He abundances 
(derived from the width of the HB) and luminosity is very good, 
there are a few clear exceptions, which are massive clusters 
with very short HBs (and hence small values of 6Y), the clear- 
est examples being NGC 104 (47 Tuc), and NGC 5024, and 
to a lesser extent NGC 5272 (M 3). We emphasize that the 
deviation of these clusters from the trend shown by the vast 
majority of the clusters is much larger than observational er- 
rors. It is notable that two of the most well studied clusters 
(47 Tuc and M 3), often considered as templates of metal- 
rich and metal-intermediate GCs for several studies of the HB, 
are found to be the exceptions rather than the rule. We recall 
that both NGC 104 and NGC 5272 have rather small values of 
IQR[Na/0], much less than observed in other clusters of sim- 
ilar luminosity. This again emphasizes the connection existing 
between the extension of the HB and the Na-0 anticorrelation, 
hence supporting the explanation we suggest that the width of 
the HB is related to spreads in the He abundances. 
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We finally recall that D'Ercole et al. (2008) and Carretta 
et al. (2009d) presented a scenario for the formation of GCs 
in which the correlation between cluster luminosity and exten- 
sion of its HB may occur naturally. According to that scenario, 
the formation of a GC is the final act of a series of events 
that begins from the violent onset of star formation in a very 
large original cloud. Star formation in this cloud continues un- 
til the kinetic energy injected by SNe and massive star winds 
causes dissipation of the remaining gas. In this scenario, this 
early population is not compact enough for the formation of a 
stable cluster, and should dissipate after the violent relaxation 
produced by the loss of the remaining gas and the ejecta of 
the most massive stars. However, once this very violent ini- 
tial phase terminates, a much quieter situation follows, where 
there is no longer a strong injection of energy in the ISM; the 
slow winds from stars of intermediate mass may then lead to 
the onset of a cooling flow, and a formation of a very dense 
and kinematically cold cloud, from which the present GC may 
form, mainly composed by second generation stars (the I and E 
populations of Carretta et al. 2009a). Part of the primordial stel- 
lar population remains trapped in the current GCs, and forms 
what Carretta et al. (2009a) called the P-population. In this sce- 
nario, it is quite natural to expect that the onset of star forma- 
tion within the cooling flow should occur earher in more mas- 
sive clusters; it will also probably stop earlier, because of the 
kinetic energy injected by the massive stars of this second gen- 
eration. Therefore, we may expect that the delay of the second 
generation, hence the typical polluter mass, is roughly deter- 
mined by the mass of the GC itself. A few exceptions may be 
easily accommodated in this scenario as objects characterised 
by a prolonged phase of formation of the primordial popula- 
tion, e.g., caused by the presence of a very extended region of 
star formation or companion clusters, as often observed among 
LMC populous clusters. 

This scenario is clearly constructed on the hypothesis of 
massive AGB stars being polluters. Were rotating massive stars 
the polluters, some important modifications would be required, 
because these stars lose most of their mass in the epoch of 
core-collapse SNe explosions, and it would be very difficult 
to produce the required cooling flow. Reproducing the trends 
considered in this paper would clearly be a challenge for this 
alternative scenario. 

8. The M 3-M 13 pair 

Since the paper by Van den Bergh (1967), the pair M 3-M 13 
(=NGC 5272 and NGC 6205) has played a very special role 
in the second parameter issue. This pair has been repeatedly 
studied to search for evidence of He-abundance variations (see 
Caloi & D'Antona 2005, and Catelan et al. 2009). They are 
bright, massive and amongst the easiest GCs observable from 
the northern hemisphere. M 3, with its large population of RR 
Lyrae, has often been considered the archetype of Oo i GCs, 
while M 13 is considered to play the same role amongst the 
clusters with very blue HB. It is interesting to revisit them 
within the scenario outlined in this paper. 

According to the basic parameters defined in this series of 
papers, M 3 and M 13 might not be so close twins, after aU. 



While the metal abundance is quite similar ([Fe/H]=-1.50 and 
-1.58), the ages might be different. Without considering the cor- 
rection for the different He abundance, the relative ages we 
considered are 0.88 and 1.02 (see Table 1). After this correc- 
tion is taken into account, this difference becomes even a bit 
larger (0.875 and 1.033). It is then possible that M 3 is younger 
then M 13, a result found in some careful determinations (see 
Rosenberg et al. 1999, and De Angeli et al. 2005), but not in 
ofliers (Marm-Franch et al. 2009; Dotter et al. 2010). This by 
itself might explain why the median colour of the HB of M 13 
is much bluer than that of M 3, without any need for a very 
large difference in Ymed (we actually obtained values of 0.249 
for M 3 and 0.260 for M 13). 

However, it has for a long time been known that these two 
GCs also have a large difference in the extension of their HBs. 
According to our analysis, the spread in mass along the HB of 
M 3 is small. Assuming a constant helium abundance, we ob- 
tain minimum, median and maximum masses of 0.624+0.003, 
0.670±0.001, and 0.716 Mq, respectively. Apart from a gen- 
eral offset of about 0.02 M© (which resembles the overall off- 
set of masses derived from HST and ground-based data in our 
analysis), these values are quite similar to those obtained in 
the analysis of Valcarce & Catelan (2008), although it is possi- 
ble that some very blue/faint stars considered by these authors 
were actually not studied by the Rosenberg et al. (2000a) pho- 
tometry. As discussed by several authors, the distribution of 
masses along the HB is probably not Gaussian (see Castellani 
et al. 2005 and Valcarce & Catelan 2008). If the interpreta- 
tion that we give for the mass spread along the HB in terms 
of variations in the He content is correct, we would correct 
the offset between median and minimum mass to 0.030 M©, 
and the corresponding spread in helium would be 5Y=0.023. 
Catelan et al. (2009) argued that such a small spread in He in 
M 3 would be detectable from Stromgren photometry and spec- 
troscopy of the blue HB stars, and concluded that a uniform He 
abundance more closely matches the observations for this clus- 
ter. This might suggest that small spreads in He such as the one 
we derived for M 3 might be an artifact of a non-Gaussian ran- 
dom mass-loss distribution. However, since the bluest HB stars 
(TefF > 1 1 , 500 K, consisting of some 5% of the whole HB pop- 
ulation of M 3) could not be analyzed by Catelan et al., it is not 
entirely clear that there is no real spread in helium within M 3. 

On the other hand, M 13 has a very wide range of both 
colours and masses along the HB. The spread in mass and Y 
abundances we obtained using Mmed are 5M=0.082 Mq and 
(5Y=0.065. Had we instead used Mmax, scaled down by 0.02 Mq 
to take into account the random term in mass lost, we would 
have got 5M=0.107 Mq and (5Y=0.083. We expect then a de- 
tectable spread in colours for the MS of M 13 [6iV - I) = 
0.033 mag and possibly as much as S(V - I) = 0.043 mag, 
if A/max were adopted, assuming that the minimum Y value is 
smaller than Y^ed- This colour spread could possibly be de- 
tected in the future by a careful analysis of HST photometry. 

As noticed by several authors, this huge difference in the 
width of the HB between M 3 and M 13 is quite unexpected, 
since M 3 is a luminous cluster (My = -8.93), even more lu- 
minous than M 13 (My = -8.70). As the discussion in Sect. 
6.4 shows, we should indeed expect a more extended HB for a 
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luminous cluster like M 3. As a matter of fact, M 3 (and more- 
over 47 Tuc and NGC 5024) appears peculiar with respect to 
the other GCs, while M 13 appears to be similar to the ma- 
jority of GCs in this respect. Within the scenario we described 
in the previous section (see also Carretta et al. 2009d), we are 
tempted to invoke a delayed cooling flow to explain the pecu- 
liarity of M 3 and 47 Tuc. This might explain other differences 
between M 3 and 47 Tuc on the one hand, and the other mas- 
sive GCs on the other, such as the small values of IQR[Na/0]. 
However, we admit that without an y f urther information, this 
hypothesis might look rather ad hoc [3 Thus, ironically, while 
the present approach seems to explain most of the observations 
relate to the second parameter effect, its application to the most 
classical couple M 3-M 13 still encounters difliculties. 

9. Conclusions and suggestions for furtlier work 

As we have seen, our reanalysis of public extensive photomet- 
ric databases of GCs demonstrates that age is the main sec- 
ond parameter affecting the HB morphology. This hypothesis 
is able to explain quite well most of the observables related to 
median HB stars of GCs, when coupled with a simple mass-loss 
law that is a linear function of [Fe/H]. Among the many observ- 
ables that are successfully explained, we note the Oosterhoff 
dichotomy that we attribute to the peculiar age-metallicity dis- 
tribution of Galactic GCs. Oo n clusters are mostly old, while 
Oo I are predominantly young, although young Oo ii and old 
Oo I GCs exist. However, at least a third parameter is required 
(and possibly even others) to fully explain the median colours 
of HBs (in particular those with very blue HBs) as well as their 
extension. There are various reasons to identify this third pa- 
rameter with variations in the He content. These include the 
variation in the scatter with metallicity, some correlation with 
the R-parameter, and the clear links with chemical anomaUes 
observed in GCs. This result is strongly indicative of a possi- 
ble link between the colours of the stars on the HB and their 
original composition, in a multiple generation scenario for the 
formation and early evolution of GCs. Self-pollution in GCs is 
possibly responsible for a large variety of the second parameter 
features, and may be in part described using the Na-O anticor- 
relation, although some modulation according to cluster lumi- 
nosity is required[3 

A combination of age and He variation therefore appears to 
be an explanation of the long-standing problem of the second 

After the first draft of this paper was written, Lane et al. (2010) 
suggested that 47 Tuc might have formed by the merging of a bi- 
nary cluster. This suggestion comes from a completely different line of 
thought, being based on the internal kinematics. As mentioned in the 
previous section, it is possible to speculate on a possible connection 
between such an origin of 47 Tuc and its anomalously short HB. 

The scenario we propose should of course not only explain the 
Milky Way GCs, but e.g., those of Fornax (Buonanno et al. 1998). This 
case is quite puzzling, since clusters 1, 3, and 5 are nearly coeval, and 
have similar metallicity, and still have very different HB's. We did not 
quantify these variations in terms of mean colours and magnitudes as 
done for the Milky Way GCs considered in the present paper, hence we 
cannot provide any quantitative analysis. We only note that the ranking 
of HBR ratios for the three coeval clusters of similar metallicity (1, 
3, and 5: -0.2, 0.50, and 0.44) is the same as the ranking in absolute 



parameter, although we do not exclude additional parameters 
such as the CNO abundances or the presence of binaries (e.g., 
of the progeny of blue stragglers) possibly playing some role. 
However, this issue is still far from being completely settled. 
We need to make some progress in developing models, and a 
number of observational tests. A short list includes: 

- Understanding the nature of the polluters. This requires ad- 
vances in the modeling of AGB stars and rotating mas- 
sive stars. Furthermore, detailed spectroscopic data for stars 
in massive and very young clusters, such as RSGCl and 
RSGC2 (Davies et al. 2007), or intermediate age clusters 
in the LMC, where multiple MS TO's have been observed 
(Mackey et al. 2008; Milone et al. 2009), may provide a 
crucial test of this scenario. 

- The present discussion has focused on He, Na and O, but 
the abundances of other elements may also play an impor- 
tant role. For instance, Al might be a proxy for He that is 
more reliable than Na. Unfortunately, our data for Al are 
not as extensive as those for Na and O, but the relation be- 
tween the production of He (which most likely occurred 
within MS stars), and the proton capture processes (which 
might have occurred in the same main sequence stars, if 
massive and fast rotating, or later during the AGB phase, if 
the stars were of intermediate mass) must be clarified. 

- A number of confirmations of this scenario are required. 
These include (i) direct determinations of He, Na, and O in 
HB stars, which were shown to be possible in at least some 
cases by Villanova et al. (2009); and (ii) a discussion of the 
luminosity of the RGB bump that we defer to another paper 
currently in preparation. 

- In addition, we ask: do properties of RR Lyrae variables 
agree with expectations? Are anti-correlations found where 
expected (important clusters such as M 54 and NGC 1851 
do not yet have adequate data)? Are multiple sequences ob- 
served where they are expected? We note here that while 
the connections between multiple MSs and variations in 
the He abundance is quite clear, the case of multiple sub- 
giant branches (SGB) is more ambiguous. SGB splitting 
measured using visual-red-near infrared colours (see e.g., 
Milone et al. 2008) might be due to a variation either in 
age or most likely total CNO content (see e.g., D'Antona 
et al. 2009 and Cassisi et al. 2008), or even [Fe/H] (in this 
case however some spread in the MS and RGB is also ex- 
pected). Interpretation of splitting is even more ambigu- 
ous when considering ultraviolet colours, which are sen- 
sitive to N excesses. Variations in He abundances are only 
marginally effective in these cases, because sequences dif- 
fering only in Y are very close each other on the SGB 
(D'Antona et al. 2002). Variations in total CNO content 
can be most likely attributed to the contribution of ther- 
mally pulsing AGB stars (Cassisi et al. 2008), which have a 
rather low mass and probably do not contribute much to He 
abundance variations. It is then unclear that there should be 
any correlation between the SGB splitting and large spreads 
on the HBs. In fact, NGC 2808 has a quite narrow SGB 

magnitude My (-5.32, -7.66, -6.82). A quantitative analysis is required 
to settle this issue. 
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(Piotto et al. 2007). SGB splitting has been detected us- 
ing visual-red-near infrared colours in 47 Tuc (Anderson 
2009), NGC 1851 (Milone et al. 2008), and NGC 6388 
(Moretti et al. 2009). These clusters have very different HB 
morphologies, ranging from very short (47 Tuc), to bimodal 
(NGC 1851), to very extended (NGC 6388). We obtain very 
different estimates of the He spread (0 for 47 Tuc; 0.048 
for NGC 1851, a possibly too large value compared with 
those determined by Salaris et al. 2008, and Catelan et al. 
2009b; 0.037 for NGC 6388). This lack of coiTelation sug- 
gests that the two phenomena are somewhat different, as 
expected if the mass range of the polluters changes from 
cluster-to-cluster. 

- The scenario requires a number of refinements. All analytic 
dependencies adopted throughout this discussion should be 
reviewed, and possibly replaced by a comparison with syn- 
thetic HBs. This may allow us to detect additional effects, 
e.g., a variation in total CNO abundances, not included in 
the present analysis. 

- Finally, hydrodynamical simulations of the formation and 
early evolution of massive star clusters are urgently needed. 
There are aspects of the present scenario that are necessary 
to explain observations, but should be understood more 
clearly. The most intriguing is the existence of a pool of 
gas from which second generation stars formed, which is 
composed of material processed through H-burning at high 
temperature diluted with pristine gas. How this pool of gas 
is generated, and how the stars form from it within the po- 
tential well of the young GC, while other stars of the orig- 
inal population evolve, remains unclear. Some explorative 
results were obtained by D'Ercole et al. (2008) which while 
very promising should be placed on a sounder basis. 
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